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ABSTRACT 



Hyperaccreting neutron star or magnetar disks cooled via neutrino emission 
can be a candidate of gamma-ray burst (GRB) central engines. The strong 
field > 10 15 — 10 16 G of a magnetar can play a significant role in affecting the 
disk properties and even lead to the funnel accretion process. In this paper we 
investigate the effects of strong fields on the disks around magnetars, and discuss 
implications of such accreting magnetar systems for GRBs and GRB-like events. 
We discuss quantum effects of the strong fields on the disk thermodynamics 
and microphysics due to modifications of the electron distribution and energy in 
the strong field environment, and use the magnetohydrodynamical conservation 
equations to describe the behavior of the disk flow coupled with a large scale 
field, which is generated by the star-disk interaction. If the disk field is open, the 
disk properties mainly depend on the ratio between \B^jB z \ and Q/Qk with B^ 
and B z being the azimuthal and vertical components of the disk field, Q and Qk 
being the accretion flow angular velocity and Keplerian velocity respectively. On 
the other hand, the disk properties also depend on the magnetar spin period if the 
disk field is closed. In general, stronger fields give higher disk densities, pressures, 
temperatures and neutrino luminosity. Moreover, strong fields will change the 
electron fraction and degeneracy state significantly. A magnetized disk is always 
viscously stable outside the Alfven radius, but will be thermally unstable near the 
Alfven radius where the magnetic field plays a more important role in transferring 
the angular momentum and heating the disk than the viscous stress. The funnel 
accretion process will be only important for an extremely strong field, which 
creates a magnetosphere inside the Alfven radius and truncates the plane disk. 
Because of higher temperature and more concentrated neutrino emission of a ring- 
like belt region on the magnetar surface covered by funnel accretion, the neutrino 
annihilation rate from the accreting magnetar can be much higher than that from 
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an accreting neutron star without fields. Furthermore, the neutrino annihilation 
mechanism which releases the gravitational energy of the surrounding disk and 
the magnetically-driven pulsar wind which extracts the stellar rotational energy 
from the magnetar surface can work together to generate and feed an ultra- 
relativistic jet along the stellar magnetic poles. 

Subject headings: accretion: accretion disks — black holes — gamma rays: bursts 
— magnetic fields — neutrinos — stars: neutron 

1. Introduction 

Hyperaccreting black hole systems formed in massive star collapses or compact object 
binary mergers have been considered as a candidate for gamma-ray burst (GRB) central 
engines for two decades (e.g., Eichler et al. 1989; Narayan et al. 1992; Woosley 1993). 
The accreting systems can drive ultra or mildly relativistic jets with huge energy via neu- 
trino (vis) annihilation process (Popham et al. 1999) or magnetohydrodynamical (MHD) 
mechanism, such as the Blandford-Znajek mechanism (Blandford & Znajek 1977) or the 
magnetorotational instability (MRI, Balbus & Hawley 1991), and finally produce GRB ex- 
plosions. The neutrino annihilation efficiency above hyperaccreting disks around black holes 
with the elaborate considerations of disk geometry, rotation, and general relativity effects 
have been studied by many authors (e.g., Ruffert et al. 1997, 1998; Popham et al. 1999; 
Asano & Fukuyama 2000, 2001; Miller et al. 2003; Birkl et al. 2007). Although the annihila- 
tion process can provide sufficient total energy up to 10 50 ergs s _1 above neutrino-dominated 
flows with accretion rate ~ 1M s" 1 (Gu et al. 2006; Liu et al. 2007), it is still under debate 
whether such a process can successfully generate a relativistic jet from the polar region of the 
central black hole. For example, MacFadyen et al. (2001) discussed that the central black 
hole formed in the "collapsar scenario" is more frequent to be formed by fallback process 
after a mild explosion (Type II collapsar) rather than formed promptly (Type I collapsar). 
The Type II collapsar can establish an accretion disk with accretion rate ~ 0.001 — 0.01M Q 
s -1 , which is not sufficient to produce a jet by neutrino annihilation. Numerical simulations 
showed that the MHD mechanism may be more efficient to drive an energetic magnetically- 
dominated wind and generate a GRB explosion (Tchekhovskoy et al. 2008, 2009; Nagataki 
2009), or the annihilation process and MHD winds can work together to provide the energy 
for GRB explosions (Harikae et al. 2009). 

On the other hand, both observational and theoretical evidences show that newborn 
neutron stars or magnetars rather than black holes can form in the GRB central engines 
(e.g., Dai & Lu 1998a, 1998b; Zhang & Meszaros 2001; Dai 2004; Dai et al. 2006; Mazzali 
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et al. 2006; Soderberg et al. 2006; Shibata & Taniguchi 2006; Lee & Ramirez- Ruiz 2007). 
Highly magnetized neutron stars formed by accretion-induced-collapse, convection, a — Q 
dynamo mechanism or differential rotation can extract their rotational kinetic energy up 
to ~ fO 51 — fO 52 ergs by spinning down via magnetic activities (e.g., Usov f992; Duncan 
& Thompson f992; Kluzniak & Ruderman 1998). A thermally- driven neutrino wind is 
dominated during the Kelvin-Helmholtz cooling epoch after the neutron star formation, 
lasting from a few seconds to tens of seconds based on the strengths of surface fields. After 
that, magnetically-dominated or Poynting-dominated wind becomes significant along the 
polar region (Wheeler et al. 2000; Thompson 2003; Thompson 2004; Metzger et al. 2007). 
The jets from magnetars may be accelerated to higher Lorentz factor ~ 100 — 10 3 at large 
radius several tens of seconds after core bounce, and bring its magnetic energy to kinetic 
energy, which is finally dissipated in regular GRB internal shocks (Thompson et al. 2004; 
Metzger et al. 2007; Bucciantini et al. 2007, 2009; Lyubarsky 2009a, 2009b). 

However, all the magnetized neutron star or magnetar models ignore the accretion 
process which occurs onto a protoneutron star for the first several seconds (Woosley & 
Bloom 2006). In the scenario of massive star collapse, the outgoing shock generated by a 
successful Type Ibc supernova explosion with a velocity of ~ 10 9 cm s _1 can evacuate a 
cavity around the new born compact star in the center of a supernova remnant. However, 
whether the outgoing shock in the core collapse of a massive star can dominate the accretion 
process and turn the accretion around is still unknown. If the rotational core collapse can 
lead to the formation of a neutron star or a magnetar, there is possible that the prompt 
accretion or fallback process can make a hyper accreting disk around the young formed star. 
On the other hand, recent simulations also showed the possibility of a debris disk around a 
massive neutron star as the outcome of a compact binary merger (e.g., Shibata & Taniguchi 
2006, Lee & Ramirez- Ruiz 2007). Similar to the black hole system, the transiently existing 
torus around a neutron star or magnetar with sufficient angular momentum can form a 
neutrino-cooled disk and release its gravitational binding energy via neutrino emission. 

The hyper accreting neutron star or magnetar system is also proposed as another possible 
central engine of GRBs. Zhang & Dai (2008a, 2009, hereafter papers I and II) found that, 
due to the inner surface boundary of the compact star, the disk has a denser, hotter inner 
region with higher pressure compared to its black hole counterpart. Also, the entire disk 
can cool more efficiently via neutrino emission and produce a higher neutrino annihilation 
luminosity. A heavily thermally-driven outflow from the surface of the new-born neutron 
star at early times (~ 100ms) prevent the outflow from being accelerated to a high Lorentz 
factor (Dessart et al. 2009). However, if the disk accretion rate and the neutrino emission 
luminosity from the surface boundary layer are sufficiently high, an energetic ultrarelativistic 
jet via neutrino annihilation can be actually produced after the jet breaks out of the stellar 
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envelope around the disk-neutron star system. 

Until now, none of the above works consider the effects of magnetic fields on the hy- 
peraccreting transient disks around the central compact stars. The MHD mechanisms in 
the accreting black hole systems mainly focus on energy extraction near the central black 
holes, and the proto-magnetar magnetically-driven wind model discusses wind generation 
and propagation along the polar regions. In the hyper accreting disk the spherical Alfven 
radius of the central star can be estimated as 

r A ~ 0.207Mr 2 2/7 M- 4 1/7 /^ 7 km, (1) 

where M = O.O1M_ 2 M s" 1 is the accretion rate, M = IAM 1A M Q is the central star mass, 
\i = /i 30 10 30 G cm 3 is the central magnetic flux. Figure 1 shows the Alfven radius r A as a 
function of accretion rate with various magnetic fields. If the surface field B is less than 

B < B ciit = 0.89 x lO^M^M^r^G (2) 

with r = r^glO 6 cm being the star radius, the accretion flow will continue to be confined 
in the disk plane without co-rotating with the compact object or getting funneled onto the 
magnetar poles. Papers I and II assumed the strength of stellar surface field is less than 
10 15 G, and did not consider the effects of a magnetic field in the disk. Recently, Xie et 
al. (2007, 2009) and Lei et al. (2009) studied the structure of magnetized hyperaccreting 
neutrino-cooled disks around black holes based on similar models, while the magnetic fields 
are generated in the disks up to 10 15 — 10 16 G (see also Janiuk & Yuan 2009). Xie et 
al. (2007, 2009) showed that the magnetic braking and viscosity can drive magnetically- 
dominated accretion flows rather than neutrino-dominated, and turn the disk temperature 
to be lower than that without field. Lei et al. (2009) investigated the properties of the 
NDAF with the magnetic torque acted between the central black hole and the disk. The 
neutrino annihilation luminosity can be increased by one order of magnitude for accretion 
rate ~ O.5M s -1 , and the disk becomes thermally and viscously unstable in the inner 
region. However, all of their works neglected microphysics processes and the changes of 
equations of state in strong magnetic fields, which may change the disk pressure and various 
neutrino cooling rates in the neutrino-cooled disks significantly. Furthermore, if we look 
into the neutron star disks, we should keep in mind that the structure of magnetic fields 
in disks around magnetized neutron stars or magnetars are very different from their black 
hole counterparts. The origin of a magnetic field in this case in the central compact star 
constructs the initial topology of the strong field. An interaction between the star and the 
surrounding disk makes the stellar magnetic field partially thread the accretion disk. 

Our motivation in this paper is to investigate the effects of a strong magnetic field on 
the hyperaccreting disk around a magnetar formed in the collapse of a massive star or the 
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merger of a compact object binary, and their observational phenomena related to GRBs and 
associated events. The accretion process onto the magnetized neutron star has been widely 
studied in X-ray binaries, T Tauri stars and cataclysmic variables since Ghosh & Lamb 
(1978, 1979a, 1979b), both theoretically (e.g., Konigl 1991; Ostriker and Shu 1995; Lovelace 
et al. 1995; Wang 1995; Lai 1998), and numerically (e.g., Hayashi et al. 1996; Miller & 
Stone 1997; Romanova et al. 2005, 2008; Long et al. 2007, 2008). The neutrino-cooled 
hyperaccreting disks, on the other hand, are much hotter, denser with much higher pressure 
and accretion rates compared to the normal disks. Therefore the structure and radiation 
process in these magnetized disks must be very different from the normal ones. For example, 
the size of the magnetosphere near the central star can be ignored except for a magnetic 
field > 10 16 G. The pressure relation B 2 /8n 3> P ma ttcr + pvl will not be satisfied in most 
cases. More attractively, we need to study the neutrino emission and annihilation process 
which is almost the most important properties of the hyperaccreting disks but never occur 
in the normal magnetized ones. 

This paper is organized as follows. In Section 2 we consider the quantum effects of 
a strong magnetic field in the microphysical scale on disk density, pressure and various 
neutrino cooling rates. In Section 3, we discuss the disk conservation equations coupled 
with the global fields both from the central magnetar and generated via the magnetar-disk 
interaction. Also, we discuss the field topology of the central magnetar. Combining the 
equations in Section 2 and 3, we discuss the properties of the hyperaccreting disks in the 
strong magnetic fields numerically in Section 4. The structure of the magnetized disks 
depends on the detailed field strength and configuration. In Section 5.1, we further discuss 
the disk outflows, which may provide the kinetic energy of the supernova associated with a 
GRB (MacFadyen & Woosley 1999; Kohri et al. 2005). Then in Section 5.2, we describe 
the funnel process onto the magnetar surface when the field strength is extremely high. 
The funnel process will significantly affect the neutrino annihilation efficiency Section 6 is 
a broad discussion. We consider the importance role played by Joule dissipation, nucleon 
reactions and r-process occurring in the magnetized disk, as well as various outcomes of the 
massive star collapse using a unified point of view. Moreover, we particularly focus on the 
application of magnetized hyperaccreting disks in the GRB and GRB-like events (e.g, X-ray 
flashes), and use a unified point of view to discuss the core collapse models related to GRBs. 
Conclusions are presented in Section 7. 
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2. Thermodynamics and Microphysics in Strong Magnetic Fields 

2.1. Density and Pressure in Magnetic Fields 

The distribution and energy of charged particles will change significantly if the magnetic 
field is greater than the critical value B ci = m 2 c 3 /qih, where rrii and qi are the mass and 
charge of the particle respectively. As the critical value for electrons is B ce ~ 4.4 x 10 13 G 



and B cp ~ 1.5 x 10 20 G, we have to adopt the relativistic Dirac equation for electrons in 
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the environment near a magnetar with the field B > AA x 10 G, while the protons can be 
still considered as classical. The main modification of electron distribution is that, the phase 
space factor in the absence of the magnetic field should be replaced by the summation over 
possible Landau level 

T3 / d 'P — > JLdnL / -T^dp, (3) 
J n=0 J 

where B m is the field strength, and g n L = 2 — 5 n0 with 5 n0 being the Kronecker delta function. 

The electron energies are given by (Johnson & Lippmann 1949) 

E e = ^p 2 c 2 + + 2n L eB m hc = m e c 2 \/x 2 + 1 + n L b, (4) 

where n labels the Landau level, b = 2ehB m /m 2 e c' is the dimensionless magnetic field pa- 
rameter, and x = p/m e c is the dimensionless momentum. 

The number density of electrons and positrons with the magnetic field strength B m is 

n J-OO ft C 

^ ^® POO 

\~YJ h ^l9nL j dxf e T, (5) 
~ — n ^ 



n=0 
27T 



n=0 



where 



/e=F = 1 1 (6) 

e m e c 2 ^x 2 +l+bn L /k B T^fri e _|_ ^ 

is the Fermi-Dirac function. The total pressure in the disk is contributed by five terms, i.e., 
the pressure of nucleons, electrons (including positrons), radiation, neutrino and magnetic 
field: 

P = P nuc + P rad + P e + P u + P B . (7) 

We take the approximation that P nuc and P ra d do not change with magnetic fields^, and 
P v = u u /3 is only noticeable in very opaque regions in the disk, where u v is the energy 



^^We neglect the photodisintegration process which happens far from the central star r > 400 km, thus 
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density of neutrinos. The pressures of electrons and positrons in the strong magnetic field 
are 



and the electron pressure P e = P e - +P e +- The pressure of the magnetic field is Pb = -B„/87r. 

If B — y 0, the number density formula ([5]) and ([8]) as a series of different Landau levels 
will switch back to their integral form without fields. However, the convergence performance 
of Landau level series becomes poor for weak magnetic fields. Figure 2 shows the convergence 
performance of density and pressure Landau level series, and compares the summations of 
formula ([5]) and ([8]) to the values without fields. Here we take temperature and chemical 
potential as fixed typical values in neutrino-cooled disks. We define b n as the n term in the 
Landau level series of §5§ and ([!])• We note that the value of b n drops faster for stronger 
magnetic fields. So we have to take a huge number of b n terms to calculate the summation 
for relatively weak magnetic fields. Moreover, whether the microphysics change in the strong 
fields is important also depends on the disk temperature. The difference between density or 
pressure with and without fields is more significant for lower temperature. This conclusion 
can also be applied to the hyperaccreting disks around black holes which generate the fields 
by MRI or different rotation. Later in this paper, we take the upper limit n = 10 4 for b n 
summation calculation for B > B ce , as the pressure or the equations of state in a magnetized 
disk switch back to those without fields for B < B ce and the typical temperature in the disks. 



Neutrino cooling processes (especially, the Urea process) in the environment of compact 
object magnetic fields have been studied for years (e.g., Chen et al. 1974; Dorofeev et al. 
1985; Dai et al. 1993; Yuan & Zhang 1998; Lai & Qian 1998; Roulet 1998; Baiko & Yakovlev 
1999; Yakovlev et al. 2001; Duan & Qian 2004, 2005; Luo 2005; Riquelme et al. 2005). 
Similarly, strong magnetic fields also affect the neutrino emission in a hyperaccreting disk 
around the compact object. Here we consider the problem systematically in this section. 



the nucleons are mainly composed of classical protons and neutrons without a-particles in the disk region we 
focus on. The equations of state for protons and neutrons do not change with field B < 10 20 G as mentioned 
above. 



P 
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2.2. Neutrino Cooling in Magnetic Fields 
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The total neutrino cooling rates in the vertically-integrated disk are taken as (Di Matteo et 
al. 2002; Kohri et al. 2005) 

Q - = y ( 7 / 8 >^ r4 (9) 

where the total optical depth for three types of neutrinos r n = r a jy . + T s>Ui with r a i/ . being 
the absorption depth and T a>Vi the scattering depth, both of which can be affected by strong 
magnetic fields. 

The neutrino absorption depth can be approximately given by (Popham et a. 1999) 
T a)U . = q a ,u t H I '(!<t bT a ) , where q a>v . are the absorption neutrino cooling rates for three types 
of neutrino, and H is the half thickness of the disk. The electron neutrino cooling rate q a ^ e 
in the disk can be simply taken as the summation of four terms q a ^ e = (q e N + q e -e+^u e u e + 
(jVems + QpiamsorO-ff, which are the the cooling rates due to electron-positron capture by 
nucleons, electron-positron pair annihilation into neutrinos, nucleon-nucleon bremsstrahlung, 
and plasmon decays respectively (Kohri & Mineshige 2002). On the other hand, the muon 
and tau neutrino cooling rates are the summations of pair annihilation and bremsstrahlung 
qa,u„H = q a ,v T H ~ (q e - e +^u T u T + qbmms)H. 

The neutrino cooling via electron-positron capture by nucleons (or say the Urea process 
in the disk), which is usually the most significant cooling rate among various cooling terms 
in hyperaccreting disks, can be considered as the combination of three processes: q e N = 
q P +e-^n+u e + q n +e+^ P +u e + q n ^ P +e-+D e - If we take the parameter K as 

k = g|gKj_+ ggi) (mec)6; (1Q) 
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where G\ ~ 1.436 x 10" 49 ergs cm" 3 , C v = 1/2 + 2sm 2 6 w with sm 2 6 w = 0.23, and g A = 
— 1.23. Using the modification rule in strong fields ([3]), we can obtain the three terms of 
neutrino cooling rates in strong magnetic fields as 

q p+e -^ n+VB = b— > g nL n p / =?/«-» ( n ) 

4 ^ Jmax{q,VT+^Zb} V £ ~l~ Tl L b 

h K^ [°° eje + qfde 

q n +e+^p+v e = 0— 2_.9nLn n I - =r/e+> { 12 > 



qn^-p+e +u e 

n=0 



4± 9 „ L n n r 7 t7> 3& (1 -/.-), (13) 
4 ±^ JVT+^b V? 2 -l-n L b 
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where e = E Ue /m e c 2 is the dimensionless energy of electrons, and q = (m n — m p )/m e ~ 2.531. 
Here we adopt the approximation adopted in Shapiro & Teukolsky (1983) that the total 
reaction on nucleons is directly proportional to the density number of nucleons, i.e., protons 
and neutrons. Figure 3 shows the convergence performance of Landau level series of two 
dominated cooling rates (jTTj) and (|T2j) with the typical temperature and chemical potential 
in neutrino-cooled hyperaccreting disks. Similar to the density and pressure in a strong 
magnetic field, the convergence properties of neutrino cooling series also go poorer for a 
weaker field. Moreover, the total Urea neutrino cooling rate becomes lower for a higher 
magnetic field. Although the difference is not significant for T = 5 x 10 10 K, it can be 
greater for a lower temperature, as shown in Figure 2. We will discuss this effect later in 
Sections 4 and 6. 

The electron-positron pair annihilation rate without magnetic fields can be calculated 
as (Burrows & Thompson 2004) 

q e -e+->v & u e - qu e p e (B=o) ^ 2.558 x l0 33 T^f(r] e ) ergs cm" 3 s" 1 , (14) 
q e - e +->v^ = qe-e+^u T u T ^ 1-093 x lO 33 !^/^) ergs cm" 3 s" 1 , (15) 
where the function f(i] e ) is 

fiVe) = 2F 4 (0)F 3 (0) (16) 

with the function F n (±r) e ) being 

FJ±rj e ) = / -dx. (17) 

In the case of strong magnetic fields, we can approximately take the annihilation rate as 

B m ehm 2 e c 5 
Q-vm — qu t Di(B=o) x 287r 2 T 4 A; 4 ' 



4J,4 ' 
B 



when the field strength B m satisfies B m hec ^> T 2 k 2 B (e.g., Kaminker et al. 1992; Yakovlev 
et al. 2001). 

In most cases the cooling processes via bremsstrahlung and plasmon decay are much 
less significant than e~e + capture and annihilation, furthermore the unchanged distribution 
and energy of neutrons and photons. Therefore we still use the formulae as in Kohri et al. 
(2005) for these two processes. 

Besides the absorption depth, three types of optical depth for neutrinos through scat- 
tering off nucleons and electrons are given by 

= (a^pYp + <r vtn Y n + a Vie Y e )pH/m B , (19) 
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where a UiP , a Uifl and a UiC are the cross sections of scattering on protons, neutrons and elec- 
trons, p is the density of the disk. We take the cross sections of a UiP and o v . n in strong 
magnetic fields as 

a np ~ a Uip(B=0) = f [j^j 2 [(C V - l) 2 + 5g\(C A - l) 2 ], (20) 

^-^o) = T ^J (—6— J' (21) 

which are still classical (Tubbs & Schramm 1975; Burrows & Thompson 2004). Here a c± 
1.705 x 10~ 44 cm 2 is the neutrino cross section coefficient. On the other hand, neutrino- 
electron scattering process in a strong magnetic field has been discussed by many authors 
using various approaches (e.g., Bezchastnov & Haensel 1996; Kuznetsov & Mikheev 1997, 
1999; Hardy & Thoma 2000; Mikheev & Narynskaya 2000). Although the analytic results 
are quite complicated and different based on various approximated treatments, all of these 
works show that the cross section a Uie is proportional to the field strength oc eB m . For 
simplicity, we use the factor 2eB m hc/ '(ksT) 2 to compare the importance of electron energies 
in strong fields and the thermal energy, and consider the cross section a u . e in strong magnetic 
fields as 

o Vi e ^ <r Vi e(B=o) x [1 + 2eB m flc/ (k B T) 2 ] (22) 

with 

- x O fe) + t) + c ^ + < c " - c ^ < 23 » 

3. Conservation Equations in Magnetized Disks and Field Topology 

3.1. Conservation Equations 

In Section 2 we discussed the microphysics and thermodynamics equations in strong 
magnetic fields, and compare them with the equations without fields. Now we consider 
the basic conservation magnetohydrodynamical (MHD) equations for a vertically-integrated 
steady-state magnetized disk. We modify the hydrodynamical equations without fields by 
adding the coupling of a large-scale magnetic field. In this case we consider the accretion flow 
is still constrained in the disk plane. The disk with an outflow structure will be discussed 
in Section 5.1, and the funnel accretion process in an extremely strong field in the stellar 
magnetosphere will be discussed in Section 5.2. First of all, the mass continuity equation in 
a vertically-integrated disk does not change 

M = -27rr£t; r , (24) 
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where S = 2pH is the disk surface density and v r is the radial velocity of the accretion flow. 

The vertically-integrating angular momentum conservation equation with fields reads 
(e.g., Lovelace et al. 1987; Shadmehri & Khajenabi 2005) 

M— = -2tt— r 3 Zv— - r 2 B z B^\ z=H + { r 2 B r B^)\ z=H - — (Hr 2 {B r B^)), (25) 
ar ar \ ar J ar ar 

where (...) = f^ H dz(...)/2h, v = ac s H is the kinematic viscosity coefficient in the disk 
with a being the viscosity parameter. Here we still use the a-prescriptionB. The magnetic 
field can play an important role in transferring the angular momentum. If the radial field B r 
can be neglected (as will be shown in Section 3.2), we can integrate the angular momentum 
equation as 



M „ MN B , . 

3^ = " E -3^ (26 > 



where 



OG 



MN B (r) = - / r l B z B^\ z=H dr (27) 



is the integrated magnetic torque, / = 1 — lo/(r 2 Q) with Iq being the specific angular mo- 
mentum constant. In the standard assumption that the viscous torque is zero at the inner 
boundary of the disk r* + b with r* being the stellar radius and () < r„ we can take 
f = 1- ^fr~Jr (Frank et al. 2002). 

The local energy conservation equation in a magnetized disk is 

Qtm + Qjoule = Qadv + Qu > ( 28 ) 

where Q^ is and Q~ Ay keep the same as in the case of B m = ( (IA12[) and (IA13[) in Appendix 
A), and the neutrino cooling rate Q~ in strong magnetic fields have been discussed in Section 
2.2. In particular, the Joule dissipation (or say the Joule heating) term Q^ oule is 



AttH H 
<2 Lie = ^(J 2 ) = -r ?t ((Vx5 m ) 2 ) 

_ Hth f 2B r \ z=H dB z _ H+ fdBzY + ]_ 
An 1 H dr z ~ \ dr ) r 2 



d(r 



dr 



(29) 



2 Actually, the half thickness of the disk H and the kinematic viscosity v = ac s H will change depending 
on the field structure. For example, the half thickness H will increase for large B§ if B r <g; B$,B Z (Zhang 
& Dai 2008b). Also, Shakura & Sunyaev's scenario (1973) v — ac 2 s /Q,K should be modified as v — a'c^/flK 
with a' being the function of magnetic fields. In this paper we use the following treatment: calculate the 
magnetic pressure Pb into the total pressure, since the isothermal sound speed c s oc P 1 / 2 , c s and the disk 
half thickness can increase for strong fields. We still adopt the classical relation H = c s /£Ik and v = ac 2 /£!/< 
for simplicity. 
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where the magnetic diffusivity parameter r] t is adopted as rj t ~ v ~ ac s H (Lovelace et 
al. 1995; Bisnovatyi-Kogan & Lovelace 1998). In Section 6.1, we will discuss the Joule 
dissipation term in more details. 

Moreover, we need to add the charge conservation equation and the chemical equilibrium 
equation, which are the same in paper II (or see Appendix A (1A7I) and (1A14j) ). 



3.2. Field Topology 

The magnetic field of a neutron star or magnetar threading the disk has the vertical 
component 

B z (r) = B (^y. (30) 

The actual magnetic field in the vertical component near a magnetar may be a mix of 
monopole {n = 2) and dipole (n = 3) fields (Thompson et al. 2004). The dipole field drops 
faster than the monopole field along the disk radius, thus the disk properties will be quite 
different with the other cases. We focus on the dipolar form n = 3 in this paper. 

The differential rotation between the disk and the star will generate a toroidal field 
(Wang et al. 1995). Following Lai (1998), we consider two possible field configurations in a 
disk, i.e., the stellar magnetic field threads the accretion disk in a closed configuration (i.e., 
the classical configuration as in Ghosh & Lamb 1979a), or the magnetic field becomes open 
(e.g., Lovelace et al. 1995). The generated azimuthal component of the magnetic field in a 
steady-state disk is 

B^\ z =h = -0B g (31) 
for an open magnetic field in the disk, and 

B^\ Z=H = -(3 (-^j B z (32) 

for a closed magnetic field, where is a dimensionless parameter, and Qk are the stellar 
surface angular velocity and Kepler velocity respectively. The radial component B r generated 
in the disk reads 

B r \ z =H oc (f^j B Z \ Z=H (33) 

In most cases, we consider v r -C v k = Qi<r and B r \ z= n ^ B z \ z= h- Thus the field strength 
B m in the disk can be obtained as B m = B z + B^ + B% m J B z + £>|. 

According to Lovelace et al. (1995), if the angular velocities of the star and disk differ 
substantially, the magnetic field lines threading the star and the disk undergo a rapid inflation 
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so that the field becomes open. As a result, the outer part of the disk will maintain an open 
field configuration, while the inner disk in the magnetosphere still has closed lines. Since 
in the hyperaccreting case the magnetosphere is very small except for an extremely strong 
center field, in which the disk plane will be disrupted and truncated in the magnetosphere 
(see discussion in Section 5), the magnetized disk where still has a plane geometry would 
be more favorable to maintain an open field. However, we still consider the two cases as 
a broader consideration. If the magnetic field in the disk is open, the magnetic torque in 
equation (1251) becomes 

MN B = (3B 2 Z ^—^, (34) 
and the integrated angular momentum equation fl26|) can be written as 

— — / + -'-VGMB 2 - = 2a— -r 3 , 35 

3n 3ir s 2n — 3 p 1 ' 2 

where we take the ratio of the angular velocity of the disk flow and the Keplerian velocity 
as a constant s = VL/VLk- For the case of a closed disk field, on the other hand, using (1321) 
we have 

MN B = PsBl^ - ^.BKGM)-^^^ (36) 



and the angular momentum equation (1261) reads 



GMM \/GA£ — — ^2 5/2 _ = 2a ^ A (37 ) 

3tt j 3n 2n-3 z 3s P r 2n-9/2 p 1 / 2 ' V ; 

where P r is the spin period of the magnetar. 

Moreover, we obtain the energy conservation equation ( 128]) in disk with the detailed 
field topology structure as 

3GMM „ Hri t r „ . x „ 9l „ 
for an open magnetic field, and 

2 



3GMM | Hth j :: - 2 
8nr 3 Aixr 2 

for a closed field. 



s(n - 1} -^r-2 



P 2 } B 2 Z = Q; dv + Q- (39) 



We make a brief summary in the end of this section. We consider the quantum effects 
of strong fields on the disk thermodynamics and microphysics in Section 2.1, and the MHD 
conservation equation with coupling of a large scale disk field in Section 2.2. The strong 
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magnetic field can change the electron distribution and energy significantly, and also change 
the cross sections of the various neutrino reactions. The large scale field can play an impor- 
tant role in transferring angular momentum besides the viscous stress, and heat the disk as 
well. In order to see the differences of disks with and without field (or more exactly to say, 
with fields < 10 14 G in the hyperaccreting disk) clearly, we list the basic equations without 
fields in Append A for a comparison. 



4. Numerical Solutions of Magnetized Neutrino-Cooled Disks 

We adopt the a turbulent viscosity model of Shakura & Sunyaev (1973) and fix a = 0.1 
in all of our calculations. First we want to study the direct quantum effects of strong 
magnetic fields on the disk, and the macrophysical field coupling in the MHD conservation 
equations independently. Figure 4 and Figure 5 give the independent results. In Figure 4, 
we keep the conservation equations as the normal hydrodynamical equations without fields 
(i.e., equations ( ]A8[) . ( lAlOj) to (IA14|) in the Appendix), and only adopt the equations of state 
fll]) to (1231) in Section 2 to see the quantum effects of a field. On the other hand, in Figure 
5, we keep the microphysical and thermodynamical equations as in the case of no field (i.e., 
equations flAlj) to (IA7I) in Appendix A and the equations in Section 2 without fields), but 
only discuss the field effect on angular momentum transfer (|2"6T) and energy heating (I28p . 
As discussed by Duncan & Thompson (1992), in principle the magnetic fields as strong as 
3 x 10 17 (P/ms) _1 G can be generated in magnetars as the differential rotation smoothed by 
growing magnetic stresses. Kluzniak & Ruderman (1998, see also Ruderman et al. 2000) 
argued that the energy stored in the differentially-rotating magnetars can be extracted by 
the process of wounding up and amplification of toroidal magnetic fields inside the star up to 
~ 10 17 G, then the ultra strong field will be pushed to and through the surface by buoyancy 
force. Based on these considerations, we adopt the magnetar surface field up to ~ 10 17 
G, and discuss the field strength from 10 14 G to 10 17 G with a dipolar extension in our 
calculations. The steady state disk can be considered as the transient hyperaccreting disk 
with a fixed time, or say, with a fixed accretion rate. On the contrary, for a time-dependent 
disk model (e.g., Janiuk et al. 2004, 2007 without fields), we have to consider the evolution 
of the magnetar surface field, which is beyond the purpose of this paper. Moreover, we do 
not consider other quantum processes in strong fields such as the generation of strong electric 
field and electron/positron pairs as E — > e~ + e + + E and j + B — > e~ +e + + B (Usov 1992). 
These processes are important in the pulsar magnetosphere, but plays a secondary role in 
the disk region outside the Alfven radius. In Figure 4, a stronger magnetic field makes the 
disk be thinner, cooler with lower temperature and neutrino luminosity. On the contrary, 
Figure 5 shows that the disk will be denser with higher pressure and neutrino luminosity, 
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and becomes hotter at least in the inner disk region near the magnetar, because the stronger 
field plays a more significant role in transporting the disk angular momentum and heating 
disk by Joule dissipation^. As a result, the quantum effects in strong fields and the field 
coupling in MHD equations play two opposite roles in changing the disk properties, one to 
decease pressure, density and luminosity with increasing field strength, and the other to 
increase them. These two competitive factors work together to establish the actual structure 
of the disk. 

In Figure 6 we consider the quantum and coupling effects together. This figure is for 
the disk structure and neutrino luminosity with accretion rate M = O.1M s _1 and an open 
magnetic field configuration in the disk. From equation (135]) we know that the ratio between 
the two parameters (3/s determines the angular momentum transfer by magnetic field. As 
the result of Figure 6, a higher ratio of (3/s leads to a higher density and pressure in the 
entire disk, higher temperature and neutrino luminosity in the inner region of the disk, as 
well as lower electron fraction at the disk radius ~ 30 km. Moreover, electrons in the disk 
becomes more degeneracy around ~ 20 — 40 km for higher (3/s, but quickly changes to be 
nondegeneracy at the inner edge of the disk toward the magnetar surface. 

Figure 7 gives the magnetized disk for various surface vertical filed Bo and open magnetic 
field configuration in the disk. Stronger fields give higher density, pressure, temperature and 
neutrino luminosity. This figure shows that the magnetic field coupling in the MHD equations 
of disk flow is more important than the quantum effects of strong fields in a microphysical 
scale. In the case of B = 10 1 ' G, the disk has no steady state solution in the inner region 
of the disk where the flow should co-rotate with the central star and further be channeled 
onto the magnetic pole along the field lines. We will discuss the funneled flows later. On the 
other hand, for B Q = 10 14 G, we stop calculation where the disk field is too weak for us to 
consider the quantum effect in strong fields as discussed in Section 2. In this case, the disk 
properties will be very similar to the case of B = 10 14 G without quantum effect even in the 
case of Bq = 0, which has been showed in Figure 5. The electron fraction decreases inward 
for B < 10 16 G, which is similar to the results without fields (e.g., paper I, Fig. 7), but a 
field ~ 10 17 G will change the electron fraction distribution significantly. The distribution 
of electron potential r] e along the radius also changes a lot in strong fields. Larger peak with 
more degeneracy electron state is obtained by stronger fields, while the change of rj e becomes 
less obvious for relatively weaker fields. 



3 In this section, we extend our calculation to the inner disk region inside the Alfven radius This 
extension is good for us to see the disk properties of magnetized disks, although we need to use another 
treatment to discuss the disk structure in the magnctosphere with an extreme strong field. 
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Above we focus on the disk accretion rate of ~ O.1M s -1 , which is typical for a 
magnetar hyperaccreting disk. Moreover generally, we can study the magnetized disks with 
different accretion rates. We discuss the disk properties with open disk field and different 
accretion rate M = 0.02,0.1 and 1 M s" 1 in Figure 8. The values of density, pressure 
temperature and neutrino luminosity become higher for higher accretion rate. This result is 
similar to the case without field. Also, lower electron fraction and more degeneracy electron 
state can be obtained for higher accretion rate, except for the inner region where the Joule 
dissipation becomes more significant than the viscous heating. Therefore, the accretion rate 
plays a very similar role in changing the disk properties to that without fields. However, we 
should point out that, there is no steady-state solution for M = O.O2M s" 1 in the inner 
disk region Bq = 10 16 G. This is reasonable and similar to the case of M = O.1M s _1 and 
Bq = 10 17 G, as such a region is already inside the Alfven radius and is almost belong to the 
magnetosphere. The extension calculation in this section cannot reach this region. 

Figure 9 shows the disk structure for a closed magnetic field in the disk. We discuss it 
besides the open disk field configuration for completeness. In this case, the disk properties 
not only depend on the disk field structure and angular velocity (s), but also depend 
on the spin period of the central magnetar. A normal pulsar has a spin period around 1 
s ~ 100 s, while the new born pulsar as a candidate for GRB central engines usually has 
a much shorter period with a timescale about tens of milliseconds or even less. Here we 
consider the magnetar with a period P r = 5, 10, 100 ms. A rapidly rotating magnetar can 
act an extra torque on the disk. However, we only consider the magnetic field torques from 
the magnetar for simplicity. A shorter period of the central magnetar decreases the disk 
density, pressure and electronic chemical potential slightly, but increases the temperature 
at 20 — 40 km and the electron fraction at r > 20 km (note that there is a peak of r/ e ) at 
r = 20 km for B = 10 16 G. The effects of period change on the disk are only obvious for the 
surface vertical field Bq > 10 16 G in Figure 9. However, since the neutrino luminosity does 
not change significantly even for Bq ~ 10 16 G, we cannot expect any obvious observational 
events from the disk plane related to the effect of period on the disks directly. 

Figure 10 shows the M — £ curves for a given disk radius r = 40 and 80 km. In order 
to see the stable performance clearly, we extend our calculation to M = 1OM s -1 , although 
such an accretion rate is impossible for an accreting magnetar system. The condition for 
viscous stability is dM/dY, > 0. Therefore the disk is viscously stable for B < 10 16 G, 
but becomes unstable in the inner region r < 40 km for Bq ~ 10 17 G with both open and 
closed disk fields until M > 1M s -1 . However, because the Alfven radius r A 3> 40 km 
for M < 1M s _1 and Bq ~ 10 17 G, such an entire instable region actually interacts in the 
stellar magnetosphere. Therefore we can conclude that the accretion flow in the disk plane 
can always be viscously stable. This conclusion is consistent with that in the case of no fields 
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(Narayan et al. 2001; Di Matteo et al. 2002; Kawanaka & Mineshige 2007), but is somewhat 
different to that in Lei et al. (2009), who showed that the disk is always unstable in the 
inner disk region. The main difference is caused by different magnetic structures between 
neutron star disks and their black hole counterparts. 

The black hole hyperaccreting disk without fields is thermally stable, at least in the 
region where the nucleon gas pressure dominates over the total pressure (Narayan et al. 
2001; Di Matteo et al. 2002, Kawanaka & Mineshige 2007). However, a magnetic field in 
the disk can make the disk be thermally unstable. We roughly estimate it with an analytic 
method. The general condition for thermal stability can be taken as (Narayan et al. 2001) 
(rflnQ + /(ilnif)|s < (d\nQ~ /d\iaH)\s. The ratio of viscous and Joule heating is 



As a result we have 



^^| S< ^| E<1 + ^| S (41) 

i.e., the value of (d\iaQ + / 'd\aH)\^ increases by unity if <3j oul dominates over Q^ is in the mag- 
netized disk. Now we calculate the ratio (dlnQ^ is /dlnH)\^ in the magnetized disk. Taking 
the open field configuration as an example. The ratio of angular momentum transferred by 
magnetic torque and viscous stress is 

^ = ( c\ ~ o\ / — — — • 4-2 = V-Ws-'M-^Mzlr'rt^^Bl^. (42) 
J vis s{2n — 3) y/GMMf r 7 / 2 

Furthermore 

^gf b = 2 (l + f- ) = 2 + 24.6/3 S - 1 M- 4 1 / 2 Mr 1 1 /- 1 < 6 r 6 7 / 2 5 2 16 . (43) 

We can see that in the outer part of the disk, the ratio (j4~3]) is still ~ 2, but it can be 
significant increased in the disk inner part where the magnetic field is sufficiently high. If 
we take the radius r ~ or r 6 ~ 2MZ^ 7 M7^ 7 Bq^ 6 tI 2 / 7 , equation fj4~3l becomes 



rflnQ 



VIS 



dlnH 



2 + 2.174/3s- 1 (l - 0.707My i 7 J B " 1 2 f)- 1 . (44) 



On the other hand, for the neutrino-cooled region, the neutrino cooling rate under strong 
fields can be approximately taken as Q~ oc ST 4 oc if 8 , or (dlnQ~ / dlnH)\j] — 8. Based 
on the above calculation, whether the gas dominated disk region is thermally stable or not 
depends on the disk radius, magnetic field strength, disk angular velocity and accretion 
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rate. For typical values (3 ~ s,M_i = l,-E>o,i6 = 1, and r is comparable to r^, we have 
(dlnQ + /dlnH)\j^ ~ 10.4, which shows a thermally unstable performance. Furthermore, The 
radiation dominated region of accretion disks are expected to be thermally unstable in the 
a-model without a large-scale field (Lightman & Eardley 1974; Shakura & Sunyaev 1973; 
Narayan et al. 2001), therefore such a region in the magnetized disk around a magnetar will 
be always thermally unstable. On the other hand, Hirose, Krolik & Blaes (2009) showed the 
disk thermal stability based on the 3-dimensional radiation MHD simulations of a vertically 
stratified shearing box, in which the turbulent magnetic field is generated by magneto- 
rotational instability. The main reason of that difference is based on different microphysics 
consideration. The stress is assumed to be proportional to the radiation pressure in the a- 
prescription, while the stress fluctuations precede pressure fluctuations, or say the magnetic 
energy fluctuations drive pressure fluctuations in the MHD simulation. However, their MHD 
simulation focuses on the accretion rate near a fraction of the Eddington limit with normal 
radiation process rather than neutrino cooling. In our paper we still use the a-prescription 
with the large-scale field. Thus we have a thermal unstable conclusion in the radiation- 
dominated regions of the magnetized disks. 



5. Disk Outflows and Funnel Flows 



So far we study the structure and luminosity of the magnetized disks with an open or 
closed field configuration constrained in the disk plane. If we define the Alfven radius ta 
as the point at which the magnetic pressure Pb = equals the ram pressure of the 

accretion material pv%, and the magnetospheric radius r m as the point where the magnetic 
pressure is approximately equal to the total matter pressure Pb — -Pmattcr +pv%, the disk plane 
will be disrupted near the magnetosphere with r < as shown in equation (1) and Figure 1, 
and the disk flow will be prevented from accretion in the disk plane at r m . Then almost the 
entire disk will be tunneled along the magnetar field lines onto only a small fraction region 
of the stellar surface at r < r m . The value of r m depends on the details of disk-magnetar 
interaction. In the normal case such as accretion onto magnetized neutron stars in X-ray 
binaries, we have P ma ttcr <C pvl- Thus r m can be calculated as Pb ~ pv^. It is estimated 
that r m is close to or slightly less than the Alfven radius (Lai 1998). Let us see the situation 
in the magnetized hyperaccreting disks. Differently we have pv% < P m atter or v r < c s in the 
disk plane for the hyperaccreting case. Figure 11 gives the pressure ratio of magnetic stress 
and matter stress, and the cooling efficiency along the radius. In this Figure we confine 
the calculation in the disk plane without considering the funnel effect. We find that the 
pressure ratio keeps -Pb /-^matter < 1 in the disk plane, thus the accretion flow can hardly be 
lifted by the field pressure above the plane. However, since the accretion flow begins to co- 
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rotate with the magnetar surface inside the Alfven radius, the viscous heating and angular 
transport which are generated by differential rotation in the disk become ignorable. Also, 
the generated azimuthal component of the disk field, which is proportional to the difference 
between the disk and magnetar surface angular velocity, should also be ignorable. As a 
result, the disk begins to cool and loss its angular momentum. The matter pressure P m attcr 
drops at the radius inside the Alfven radius. In addition, as the magnetic pressure increases 
as the accretion flow moving towards the magentar surface, the ratio P-q /^matter keeps to 
increase and can be greater than unity as well. Thus the magnetosphere (with its edge at 
r m ) can actually form inside the Alfven radius, i.e., we still have the magnetosphere region 
with r* < r m < in the magnetized hyperaccreting disks. 

Furthermore, the strong field in the Alfven radius would cause a magnetically-driven 
outflow from the inner disk region; the interaction between the disk and central star is 
favorable for launching an X-type wind from the disk-magnetosphere boundary (Shu et al. 
1994), both of which also decrease the pressure of the accretion matter and increase the ratio 
-PB/-Pmatter- On the contrary, it is most likely that an thermally-driven outflow can form in 
this outer region r > r&- Generally speaking, the thermal outflow can take away angular 
momentum and energy outside r^, make the direct flow accretion rate into the magnetosphere 
and magnetar surface be low and give a possible energy source for a supernova explosion 
associated with a GRB. 

In Section 5.1, we discuss the disk outflow caused by thermal heating. In Section 5.2 we 
discuss the funnel flow process due to the co-rotation and magnetic pressure near the stellar 
magnetosphere. 

5.1. Disk Outflows 

The thermal wind is induced by viscous heating (MacFadyen & Woosley 1999), or 
neutrino heating (Metzger et al. 2008b), depending on the dominated heating and cooling 
processes in the disk. In paper II, we study the structure of the neutron star disk outflows 
around neutron stars without fields. As shown by Narayan & Yi (1994, 1995), if the adiabatic 
index 7 < 3/2 in an advection-dominated flow, the Bernoulli constant of the flow is positive 
and a thermally driven wind can be driven from the disk. In this case, the energy carried by 
the outflow is expected to feed a supernova explosion, which is associated with a GRB in the 
collapsar scenario (MacFayden & Woosley 1999; Kohri et al. 2005). This thermally-driven 
wind can also exist in the magnetized disks beyond the Alfven radius, where the matter 
pressure dominates over the magnetic pressure, and heating energy generated by viscous and 
Joule dissipation is advected inward along the disk radius. As showed in Figure 11 (right 
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panels), the neutrino cooling efficiency increases with increasing the strength of a magnetic 
field for M ~ O.1M s _1 , and the entire disk beyond the Alfven radius becomes neutrino- 
dominated for M ~ l.OM s -1 . Therefore the outflow-driven process should be significant 
below M ~ I.OMq s" 1 . In this section, we adopt the similar treatment as in Kohri et al. 
(2005) and paper II for simplicity, i.e., we only consider the viscously induced outflows and 
assume the accretion rate varies as a power law in radius for the case M < l.OM s _1 as 



U =M out (— V, (45) 

, r out I 



where M out is the initial accretion rate at the outer edge of the disk, and £ is the outflow 
index in the disk. Stronger outflows have a lower value of £. Using the outflow structure 
151) . the angular momentum equation (1261) is modified as 



1 M MN B 

= ^ - 3^n- (46) 

Here we have assumed that the outflow has the same angular velocity as the accretion flow 
outside the Alfven radius and takes away angular momentum from the disk. 

Figure 12 gives the disk structure with a viscous and Joule heating induced thermally- 
driven outflow for various outflow index £ = 0.2, 0.6, 0.9, different magnetar surface fields 
Bq = 10 16 , 10 17 G, and a fixed accretion rate M = O.5M s -1 in the radius r = 124 km (i.e., 
30R s h of a 1.4 M Q star). Similar to the results in paper II, the density, pressure and neutrino 
luminosity decrease with increasing the outflow strength. On the other hand, the ratio of 
-pB/-Pmatter and the thickness of the disk become larger for stronger outflows. Therefore a 
disk with stronger thermally-driven outflow outside the Alfven radius carries more energy 
away into the stellar envelope at a fixed radius. However, the total outflow energy strongly 
depends on the size of the region outside the Alfven radius, where generates viscous and 
Joule heating energy and induces outflows. A strong central stellar field makes the Alfven 
radius be larger and decrease the total heating and the potential outflow injection energy. 
Table 1 gives our estimate of the maximum energy carried by the outflow, which can be 
calculated as the difference between heating energy with and without filed, 
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The maximum energy injection rate decreases significantly for stronger stellar fields, which 
make the magnetically-driven wind and funnel effect be more important than the thermal 
outflow for strong fields. 

Besides the thermally- driven outflow in the disk outer region, a strong field inside the 
Alfven radius causes the accretion flow to co-rotate with the stellar field, and is possible to 
launch MHD wind along the field lines (Metzger et al. 2008a). Moreover, it is also clear that 
an X-type field configuration near the central star is favorable for launching magnetically- 
driven wind from the disk-magnetosphere boundary (Shu et al. 1994). As showed by simu- 
lations (e.g., Romanova et al. 2008), the properties of these MHD outflows depend on many 
factors such as the field structure, the accretion rate, the related viscosity and magnetic 
diffusivity, etc. The magnetically-driven winds in the disk inner region will decrease the rate 
of accretion onto the magnetar surface. Also, the magnetically-driven winds can provide 
energy to a supernova explosion together with the thermal outflow. But this issue is beyond 
the purpose of this paper. 



5.2. Funnel Flows from Disks to Magnetars 

When the magnetic pressure becomes equal or larger than the matter pressure inside 
the Alfven radius Pb > -Pmattcr, the funnel accretion process becomes significant. For the 
case of hyperaccreting disks, this process can only be important for an extremely strong 
field, which depends on the accretion rate. For simplicity, we consider the funnel process 
to take in the region between the magnetosphere and the Alfven radius, and the accretion 
in the disk plane will be truncated at the radius r m . The scenario of tunneling process is 
similar to many previous work (e.g. Lovelace et al. 1995. Fig 2, Fig 3; Koldoba et al. 2002, 
Fig 1; Frank et al. 2002, Fig 6.4). Different to the disk with cylindrical coordinates, in this 
section we use the spherical coordinates (r ; , 9, 0) with origin at the stellar center to describe 
the tunneled flow, where 9 is the angle between the axis of the disk plane and a given polar 
radius r, is the angle in the disk plane. Thus the equation for the dipole field geometry 
line is (Frank et al. 2002) 

sin 2 # . , 

ri = r d —^-, 48 
sin B 

where is the radius at which the magnetic field line passes through the disk, and G is the 
angle between the magnetic pole and disk axis. We can approximately take r rn < < ta- 

The magnetic field along the field line reads 

B p (n) = B (^j 3 (4 - 3sin 2 £) = B (^j " (a - ^sin 2 e) ^ . (49) 
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Note that the dipole field is somewhat different to the vertical component B z (r) = Bo(r*/r) 3 
in Section 3, in which r is the radius of the disk in the disk plane. The mass conversation 
along the magnetic flux line requires the poloidal part (r, 9) of the velocity v p and the density 
p satisfies 

Attd n 

K, (50) 



B p 

where K is the constant along the flux line. Near the magnetar surface, the accretion rate 
satisfies 

4irr%pJ — )=M, (51) 

where Afl is the opening solid angle of the funnel flows, v* and p* is the average velocity 
and density in this angle Afl onto the stellar surface. We use the 2tc as the total solid angle 
because the funnel flow can be accreted towards two poles. The difference between the final 
accretion rate M onto the magnetar surface and the initial accretion rate M out in the disk 
strongly depends on the strength of magnetically and thermally driven winds. Different to 
the case of Bondi accretion onto a magnetized compact star, which form an accretion column 
to cover the magnetic pole, the funnel accretion from a half-thickness disk forms a ring-like 
belt at the latitude of 9 between arcsin y/r^Jr2 and arcsin ^r*/r m to cover a part of the 
stellar surface. Recent simulations show that the geometry shape and physical properties 
of this belt (or say the "hot spot") are very complicated, depending on the detailed field 
topology and the structure of the funnel flows (Long et al. 2008, or see Romanova et al. 2008 
for a review). In this paper we take this ring-like belt as axisymmetric around the magnetic 
pole axis for simplicity. Appendix B gives our estimate of the ratio of AQ/2tt in the case of 
hyper accreting and neutrino-cooled disks. Combining (H9~j) . f l5T?j) and f lBTj) . we have 

M(27r/Afi) ( B p 



with B pif being the strength of the dipole field on the stellar surface. 

In the normal cases, the infalling funneled flow from the disk plane onto a magnetized 
compact star usually will go through a strong shock before reaching the stellar surface (e.g., 
Ferrari et al. 1985; Ryu et al. 1996; Li et al. 1996; Frank et al. 2002). However, such a shock 
is probably unlikely to develop in the hyperaccreting disks. Let us discuss it in more details. 
The Mach number in the funnel flow along the magnetic field pole can be calculated as 
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with 7 being the adiabatic index of the disk matter. Thus 
Jt 2 ~ 8tt 



\ r \ 2 / 2?r \ 2 / P p 



4vrr2 y V^P*/ V-Pmatter/ 7^pP 

= 6.30 x 10-S-Vr2 1 M 2 1 r; 6 4 J B p - 2 16 (27r/Afi) 2 (P B /P m attc r ), (54) 

where p 12 = p/10 12 g cm -3 and B p * jl6 = P p */10 16 G. For a typical hyperaccreting disk with a 
stellar field ~ 10 16 G and the accretion rate M = O.1M s -1 , we have Pb ~ -Pmatter, P ~ 10 12 
g cm -3 , and 2tc/AQ ~ 10 2 (from Appendix B), therefore the Mach number is ^# < 1, which 
shows there is no shock wave existing in the funnel flow. A strong magnetically-driven wind 
will make the density and pressure in the funnel flow decrease, but the accretion rate M 
onto the star will also decrease. Although 2n/AQ can reach to 10 3 and the factor Pe/Pmattcr 
increases for a strong field ~ 10 17 G with M = O.1M s _1 , but the decreasing factor Pi2^p*,w 
makes the Mach number hardly exceed unity. 

In the funnel channel, the flow is accelerated via the stellar gravitational force. The 
gravitational binding energy is converted to the kinematic energy of the flow, then the 
kinematic energy will be converted to the heating energy near the magnetar surface, which 
is cooled via thermal neutrino emission. The total energy release rate in the funnel process 
can be approximately taken as 

■ GMM f 1 1\ 

£funnel~— ^ [~ ~ ~ ) ^ ^5) 

and the energy equation in the magnetar surface boundary is 

loJT GMM 

8 t„ 4r, y 1 

where the parameter e is introduced to show the combined efficiency of acceleration and 
cooling. The accretion rate here is adopted as the final accretion rate onto the magnetar. 
The of the "hot spot" is S = Afir 2 . Here we consider the new born magnetar with 
a lifetime ^> 100 ms (Dessart et al. 2009), then the temperature of the "hot spot" can be 
significantly higher then the other region of the magnetar surface. The temperature in the 
"hot spot" can be calculated as 

T = 7.4 x 10 10 (27rM 1 . 4 M_ 1 e/Afi) 1/4 r ! 7 3/4 K, (57) 

where we take e as 0.5. Since the temperature is insensitive to the neutrino optical depth 
as r^ 4 , we take t v ~ 1 in our calculation. The temperature in this surface region increases 
slightly with increasing the ratio of 2?r/Afi: T = 1.3 x 10 11 K (10 MeV) for 2vr/Afi = 10, 
or T = 2.3 x 10 11 K (20 MeV) for 27r/Afi = 100. 
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The neutrino pair annihilation process i/j+Pj — > e~+e + is the most important mechanism 
to provide the energy for a relativistic jet formed from the neutron star disk with weak 
stellar magnetic field. According to paper II, the neutron star disk with a hot stellar surface 
layer could increase the neutrino annihilation luminosity by about one order of magnitude 
higher compared with the black hole disk. If the stellar field is strong up to > 10 15 G, 
as discussed in Section 3, the disk with strong magnetic field will have a higher density 
pressure and neutrino luminosity, therefore the annihilation rate above the disk plane will 
also be higher than that emitted from a disk with the same radius range without field or 
with weak field. This conclusion is similar to that in Lei et al. (2009, their Fig. 2). However, 
the strong field truncate the disk plane accreting in the inner region and will decrease the 
total annihilation luminosity from the disk. In this case, the annihilation mechanism will be 
more significant from the magnetar surface where accretes the funnel flow. Similar to the 
case without field, the total neutrino annihilation rate is contributed by three components: 
the annihilation between neutrinos both from the disk, both from the stellar surface, and 
one from the disk and the other from the stellar surface respectively. If the stellar surface is 
bright to neutrino emission, then the stellar surface will contribute to the main annihilation 
luminosity in the accreting system. This is why the neutron star accretion system will have 
a brighter annihilation luminosity than its black hole counterpart. In the funnel process, 
on the other hand, the stellar surface area where has a solid angle AQ <C 2tt and higher 
latitude 9* ~ arcsin \/r^/rA to emit thermal neutrino will be more geometrical concentrative 
and hotter than that without field, so the annihilation rate will be more efficient (Birkl et 
al. 2007). Figure 13 shows the neutrino annihilation rate as a function of height along 
the magnetic poles with different latitudes of the emitting "hot spot" area. We adopt the 
formula of neutrino annihilation as in Popham et al. (1999) and Rosswog et al. (2003). 
The different point here is that the neutrinosphere is a ring-like belt around the magnetar 
surface, but not a plane disk. The neutrino annihilation rate l vV from a higher latitude 
ring-like belt area is larger near the stellar magnetic poles, but becomes smaller far from 
the poles. That is because the shorter distance of a given point along the pole to a higher 
latitude emission area plays a key role in increasing the annihilation efficiency On the other 
hand, the less average value of the angle or larger cosdkk' between neutrino pairs from lower 
latitudes plays a more important role in increasing the neutrino annihilation efficiency at a 
large height > 20 km. Whatever, the integrated neutrino annihilation energy along the pole 
region l u pdz is much higher from the higher latitude ring-like belt than that from lower 
latitudes. Therefore, funnel accretion can accumulate more powerful annihilation luminosity 
(see the caption of Figure 13). 

In the black hole disks with strong fields, the MHD effects such as the Blandford-Znajek 
mechanism or magnetic instabilities (MRI) leads to formation of a magnetically-dominated 
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jet along the magnetic poles. Such a magnetically-dominated jet will also form from the 
magnetar in a few seconds of its formation, and the magnetar rotation energy to feed the jet. 
Therefore in a hyperaccreting magnetar system the neutrino annihilation mechanism and 
the MHD mechanism will work together to form a jet, and make the jet be more powerful 
for a GRB explosion. We will discuss this issue again in Section 6. 

Finally we list the properties of the hyperaccreting magnetar system in Table [21 taking 
the accretion rate M = O.1M s _1 near the Alfven radius or the stellar surface r* (depends 
on ta > r* or < r*) as an example. 



6. Discussions 

6.1. Joule Dissipation 

In Section 3 we discuss two competitive sets of effects to affect the disk structure and 
neutrino emission, i.e., the microphysical quantum effects and the macrophysics magnetic 
field coupling in the disk MHD equations. The quantum effects decrease the pressure, density 
and luminosity with increasing the strength of field, but magnetic coupling makes these 
quantities be higher. However, in most cases, the magnetic field coupling is more important 
than the quantum effects in a microphysical scale. We discuss this conclusion in details in 
this section. In fact, in order to see the quantum effects clearly in Figure 4, we adopt a unified 
field (i.e, without a dipolar form) in the disk, and find that this microphysical effect becomes 
more significant in the disk region at relatively larger radius. However, such a situation is 
unreal if we adopt the dipolar field from the stellar surface, thus the microphysical quantum 
effects will be less important as showed in Figure 4. On the other hand, as showed in the 
angular momentum equations (|35|) . ( 137|) and the local energy equation (|38|) . (|39|) . stronger 
magnetic fields will be more important to transfer angular momentum in the disk, as well 
as to heat the disk by Joule dissipation. Therefore, it is not difficult to understand why the 
disk with a stronger field will be hotter, denser with higher pressure and brighter neutrino 
luminosity. 

In fact, the parameter /3 which gives the relation between the vertical and azimuthal 
components of the disk magnetic field as in equation (132| i.e., the closed field configuration) 
can be calculated using the Ohm's law j = <r m (E+vx B)/c and Ampere's law V x B = Airj/c. 
We have (Lee 1999) 

o 47ra m HrVL K 

P = — ^—HrlL K = , (58) 

where a m and r\ m are the microscopic electric conductivity and magnetic diffusivity. However, 
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as argued by some authors (e.g, Bisnovatyi-Kogan & Ruzmaikin 1976; Lovelace et al. 1995; 
Bisnovatyi-Kogan & Lovelace 1997), since the accretion flow are in general turbulent, the 
microscopic r\ m should be replaced by a turbulent transport parameter rjt, which can be 
considered as being comparable to the turbulent a viscosity rj t ~ v ~ ac s H. In this case, 
equation (|58|) gives the relation (3 ~ vkoT^c^ 1 . However, this new relation has its limitations. 
As discussed by Lovelace et al. (1995), the twist \B^jB z \ can never be much larger than 
unity. Therefore, in this paper, we take /3 as a parameter rather than an obtained variable. 
Similar to Lai (1998), we consider that ~ 1 as the maximum twist of the original magnetic 
field from the central stellar surface. On the other hand, we take the magnetic diffusivity 
rjt in the local energy equation to be the order of the turbulent viscosity. As a result the 
magnetic diffusivity can be smaller than the classical diffusivity showed by Ghosh Sz Lamb 
(1979a, 1979b; see Lovelace et al. 1995 for further discussion). 

The Joule heating plays a key role in affecting the neutrino luminosity from the hy- 
peraccreting disks with small radius, because the ratio of the viscous heating and Joule 
dissipation will decrease with small radius, and the Joule dissipation quickly dominates over 
the viscosity as the main heating source in the disk. Of course all of these results only 
occur in the region outside the Alfven radius r^, although we usually give our calculations 
extended to the inner region. Figure 14 shows the neutrino cooling luminosity without Joule 
dissipation. Compared with the neutrino luminosity from the disk with Joule dissipation 
in Figures 5, 7, and 9, the main difference is that the neutrino luminosity drops quickly in 
relatively strong fields at small radius without fields, except for the case that the disk has 
relatively low accretion rate and is advection-dominated. This result is consistent with the 
above discussion that Joule dissipation plays a key role in heating the disk inner region near 
the Alfven radius. 



6.2. Application to GRBs 

As showed in Section 5.2, the funnel accretion makes the ring-like belt of "hot spot" be 
more concentrative with higher latitude than the equatorial accretion without fields or with 
moderate fields. Also, the neutrino luminosity outside the Alfven radius is also brighter for 
stronger stellar fields. As a result, the neutrino annihilation efficiency will be significantly 
increased in the hyper accreting disks around magnetars compared to the normal neutron star 
disks, and definitely more efficient than the black hole disks. However, a higher concentration 
neutrino emission region means higher temperature on the magnetar surface with more 
massive neutrino-driven winds (Qian & Woosley 1996; Dessart et al. 2009), which will 
make the jet be heavily baryon-loading. The problem is that, whether the more powerful 
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annihilation and massive mass loss rate driven by neutrino absorption can work together 
to produce a relativistic jet required for a GRB explosion? This problem is similar to that 
in the neutron star disks without fields, which has been discussed in paper II. We consider 
that a GRB-related jet forms along the stellar poles (i.e., the disk axis if there is no field), 
and only the wind materials along this axis can feed the jet and affect the bulk Lorentz 
factor of the jet. The wind ejected in a off-axis direction and that evaporated from the disk 
will not affect the jet, but only affect the potential supernova explosion associated with the 
GRB. If there is not disk around the central compact star, the neutrino-driven wind from 
the star surface can be reasonably approximated as a quasi-steady spherical outflow (Qian & 
Woosley 1996). However, the hyperaccreting accretion disk around the star can change the 
distribution of the neutrino-driven wind significantly. The entire stellar surface will be very 
hot and inject a heavily mass-loading wind along the axis after its formation in a timescale 
of ~ 100 ms, and then to cool down with the power law M wilK j oc t~ 5 / 3 . In most cases, the 
surface ring-like belt region where directly accretes the disk flow will be the hottest region 
in the disk for the time ^> 100 ms. However, most parts of the wind from this ring-like 
belt should be off-axis. The polar region that directly drives a wind along the axis is cooler 
than the equatorial ring-like belt. Even though, we use the temperature of the hot ring-like 
belt to estimate the maximum strength of the neutrino-driven wind along the poles, and 
calculate the bulk Lorentz factor of the wind along the pole with field < 10 15 G in paper 
II. A moderately or ultra relativistic jet can be produced in the hyperaccreting neutron star 
system with sufficient high disk accretion rate and bright boundary emission. 

However, the situation in a strong field will be different in some aspects. First of all, 
the neutrino-to-nucleon absorption reaction rate, which is dominated by 

u e + n ^ p + e~~ (59) 

z7 e + p — n + e + (60) 

will be significantly reduced in strong fields B > 10 16 G (Duan & Qian 2004, 2005) because 
of the quantum effects. As a result, the mass loss rate from the stellar surface which is 
proportional to the absorption rates will decrease as well. Moreover, we have to note that 
part of the newly generated e~e + particles in the neutrino absorption reaction will move 
along the closed field lines from the stellar surface rather than be injected far away from the 
central star. There are some open field lines in the surface region near the stellar magnetic 
poles (Lovelace et al. 1995), but many of them are more likely to induce winds off-axis rather 
than the neutrino-driven wind. As a result, the total mass loss rate M po iar along the magnetic 
polar region will be only a fraction of the total mass loss rate M wind . On the other hand, the 
field structure around the stellar surface will also affect the neutrino annihilation process, 
because part of the generated e~e + pairs in the reaction Ui + Ui —t e~ + e + will also move along 
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the closed field lines. Therefore, part of the neutrino annihilation energy cannot be loaded 
into the jet which propagates along the magnetic poles. The precise fraction of the total 
annihilation energy feeding the polar jet depends on the energy-momentum distribution of 
the e~e + pairs above the central star and the entire disk, as well as the ratio of e~e + plasma 
pressure to the magnetic pressure. The reliable estimate of these considerations requires 
further MHD simulations. However, since the annihilation process happens in the extended 
space above the accreting system, while the neutrino absorption reaction mainly occurs 
around the stellar surface with a sufficient high density of nucleons, the strong magnetic 
fields will play more significant role in changing the properties of the neutrino-driven wind 
rather than the annihilation process. Therefore the mass loaded in the jet will decrease more 
significantly than that of the annihilation energy. Based on these above considerations, that 
is, the neutrino annihilation efficiency is larger for a magnetized disk, such an efficiency can 
even higher for the funnel accretion process, and the mass loss along the magnetic poles 
-^wind decreases in strong fields, we can conclude that the jet along the magnetic poles can 
be accelerated to a larger bulk Lorentz factor compared with its neutron star counterparts 
without fields or with weak fields. In addition, if the magnetic pole and the disk axis do not 
overlap (0 ^ 90°), the direction of the relativistic jet and the disk rotation axis will also 
not overlap. This is an interesting topic, because if the annihilation rate from the disk is 
much less than that from the star, the jet will precess along the disk axis with a period of 
the magnetar. 

Furthermore, another significant difference between a magnetar and a non-magnetized 
neutron star in the hyperaccreting systems is that, a thermally-dominated neutrino-driven 
wind from the stellar surface may switch to be magnetically-dominated instead after seconds 
of the compact star formation (Thompson et al. 2004; Metzger et al. 2007). The condition 
of the magnetic wind can be parameterized using the term of magnetization 

a = ^K_ (61) 

M wind C 3 

where $b is the magnetic flux. The parameter a, which strongly depends on the rotational 
period of the central star, is the maximum Lorentz factor the wind can achieve if all magnetic 
energy is converted into kinetic energy, either via magnetic reconnection (Spruit et al. 2001) 
or collimation by the interaction between the wind and the the stellar envelope (Bucciantini 
et al. 2006, 2007, 2009; Tchekhovskoy et al. 2008, 2009). The critical period for <r > 1, i.e., 
the wind to be magnetically-dominated and relativistic is 

Pre < 85S 0> i 6 < 6 M^_ 5 ms. (62) 

The critical period P rc increases with increasing the stellar field, and decreases with increasing 
the mass loss rate M. The properties of the magnetically-dominated wind will be different 
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with different ranges of o for the stellar period P r < P rc . The wind is collimated along the 
magnetic poles for a < 5, but will be distributed around the direction at low latitudes and 
the equatorial plane for a > 30. For example, Bucciantini et al. (2006, 2007, 2009) showed 
that, the wind-stellar-envelope interaction may provide a viable mechanism for collimating 
the jet along the polar region. Their works are based on the consideration that a cavity 
with the radius > 10 8 cm has been evacuated by the outgoing supernova shock before the 
compact star formation. Since the wind with high a is ejected around the equatorial plane 
initially, whether or not the hyperaccreting disk with a small size 10 7 — 10 8 cm can help or 
prevent the jet collimation needs to be further studied. 

It is said that the energy of the magnetically-driven wind is extracted by the central 
star spin-down process. However, in the magnetar-disk system, if the magnetically-driven 
wind can be actually collimated along the magnetic polar region, then the stellar spin-down 
mechanism and neutrino annihilation from the "hot spot" and the entire disk can work 
together to provide energy of the polar jet and accelerate the jet to an ultra relativistic 
speed. In this general case, the bulk Lorentz factor can be estimated as 

E 

(63) 



where the energy feed rate is the summation of neutrino annihilation and the magnetic 
energy from rotational extraction, 

R = R ■ A- R 

^ J^anni \ -"-^mag 

= Lyjjfk + E Tot f mt = Ly^fvyfk) + E rot f Tot 

= ^^( £ /»*/*) + "Mr^O/rot, (64) 

with f uP , f k , f TOt being the neutrino annihilation efficiency L uP /L u , the fraction of the de- 
posited annihilation energy to provide the kinetic energy of the stellar wind, and the ratio 
of magnetic to total spin-down energy respectively. M po i ar as a fraction of M W md is the wind 
mass loss rate along the polar region (see paper II for more discussion about the mass loss 
rate). The ratio of -E anni to E m3ug is 

Eanni 2.10xl0 5 M_ 1 P r 2 r J r, 3 6 f%^V (65) 



-^mag \ J 1 



/rot 



where tj being the typical spin-down timescale. Based on the analysis in this paper, we take 
e = 0.5, /„ p ~ 0.005 - 0.01 for M_ x ~ 1, f k ~ f mt , and tj ~ 7.6S c T 2 6 (P r ./lms) 2 s@ then the 



J The annihilation efficiency f v9 for a black hole hyperaccreting disk is ~ 10 4 for M — 0.1M Q s 1 , and 
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critical value of the period P for the ratio E anQ i/ E mag being less than unity is roughly P r < 4 
ms. The critical value ~ 4 ms will even decrease if the disk accretion rate is higher with a 
higher annihilation efficiency f u „, or the spin-down timescale is longer. Keep in mind that 
this critical value is less than P rc in equation f[6"2l . which shows the wind to be magnetically- 
dominated and relativistic without neutrino annihilation. Therefore, in the extreme case of 
a millisecond magnetar, the magnetically-dominated jet extracts the rotation energy as the 
main energy source. If the stellar spin period is around tens to hundreds of milliseconds, the 
polar jet can be feeded by neutrino annihilation and magnetic energy together, and has a 
bulk Lorentz factor even higher than a. On the other hand, for the magnetar period P r > P rc 
in equation (l62l)Fl the jet formation process goes back to be thermally- driven and feeded by 
the annihilation process as discussed in paper 

nH 



6.3. Disk Nucleosynthesis and r-Process Nucleosynthesis 

Besides the GRB or GRB-like (e.g., X-ray flashes) phenomena with their associated 
supernovae, the hyperaccreting mangetar system can also generate enough 56 Ni and other 
elements from the disk for the supernova (MacFadyen & Woosley 1999; Kohri et al. 2005; 
Nagataki 2006, 2007), and produce more heavier elements in the neutrino-driven wind via 
r-process nucleosynthesis. 

Simulation works based on the collapsar scenario (MacFayden & Woosley 1999; Nagataki 
et al. 2006, 2007) showed that 56 Ni produced in a jet of the collapsar is not sufficient to 
explain the observed amount in a supernova with a duration about ~ 10 s (e.g., Mazzali et al. 



this efficiency can be higher and up to ~ 10~ 3 for the neutron star disk without fields. The value of f v t, can 
be even higher for a magnetized disk and funnel accretion as mentioned above. We take f vj? ~ 0.005 — 0.01. 
However, the f uC , value is insensitive to the final results, since the period P r oc f v & ' . The spin-down 
timescale can be approximately taken as the one in Thompson et al. (2004) and Metzger et al. (2007), i.e., 
tj oc (P r /B Q ) 2 for r A < c/O. 

5 The maximum initial spin period of newborn pulsars is still under debate. For example, Spruit & 
Phinney (1998) considered that the newborn pulsar without "kicks" can reach a long initial period ~ 100 
s. However, Heger et al. (2005) showed that the spin period of the pulsar formed a the 10 — 15M Q star is 
10 — 15 ms based on their stellar evolution code. 

6 We mention that Harikae et al. (2009) have discussed the combination effects of MHD mechanism and 
annihilation process in the accreting black hole system based on the collapsar scenario. They showed that 
the collapsar with an initial field ~ 10 10 G and angular momentum of 1.5 times the angular momentum of 
the last stable orbit can finally produce the MHD outflows and obtain the strong neutrino heating in the 
polar funnel at the same time. However, how the collapsar can generate a initial field with ~ 10 10 G remains 
an open question. 
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2006; Soderberg et al. 2008). As a result, the majority of 56 Ni with mass ~ 10 2 Mq might 
be synthesized in the accretion disk or in the disk outflows. Nagataki et al. (2007) calculated 
that the 56 Ni synthesized in the disk can reach ~ 1O~ 3 M , but it should be carried out by 
a later phase outflow or the viscous induced outflow. Also, Kohri et al. (2005) discussed 
that the recombination process of nucleons (n, p) into nuclei can happen in the outflow and 
release an energy about 8 MeV per nucleon. In the hyperaccreting magnetar scenario, on the 
other hand, since the disk becomes hotter, denser with higher pressure with a stronger field, 
the ejected 56 Ni via the disk outflow should be more than that estimated in the some region 
of the collapsar without fields. However, if the magnetar field strength is high enough and 
the funnel accretion is obvious, the disk region outside the Alfven radius r^ will be relatively 
small and ejected even less 56 Ni. Also, in the co-rotation region without viscous heating at 
radius less than r^, the disk will be cooled down and produce much less 56 Ni. 

The main properties to determine the r-process production are the asymptotic wind 
entropy S a , asymptotic electron fraction F e a and dynamical timescale td yn - In general, higher 
wind entropy, lower electron fraction and short dynamical timescale are more favorable to 
produce heavier nuclei with higher maximum A number (Meyer & Brown 1997; see also 
Thompson 2003). Although the winds are thought to be a candidate for the r-process, 
recent models without fields showed that the winds are difficult to produce robust r-process 
nucleosynthesis for the "classical" neutron star with a 1.4M mass and 10 km radius (e.g., 
Qian & Woosley 1996; Thompson et al. 2001). Thompson (2003) roughly estimated the 
r-process in the strong field environment with B > 6 x 10 14 G. Since the strong field can 
trap the wind in the neutrino heating region, the amplification of the trapping timescale 
can lead to the amplification of the entropy. This amplification may be sufficient to yield 
robust third-peak r-process nucleosynthesis. In addition, Metzger et al. (2007) showed that 
the presence of a magnetar field ~ 2 x 10 14 G and mildly rapid rotation ~ 10 ms moves 
the ratio (5 ,£l ) 3 /td y n ! which is the critical wind parameter to determine the condition for 
an r-process, an order of magnitude more favorable for third-peak r-process nucleosynthesis. 
Therefore, the magnetic fields from the central magnetar may play significant role in inducing 
a successful strong r-process nucleosynthesis. However, the problem is that, all of the works 
considering magnetic fields do not include the strong field quantum effects on the state of 
equations of the wind. In other words, they do not consider the Landau level effect and 
the neutrino absorption reaction rate modified by the strong fields. As discussed in Section 
6.2, the strong field can decrease the neutrino absorption reaction rate, and also decrease 
the mass loss rate of the neutrino-driven wind as well. Also, the strong field can affect 
the density and pressure distribution in the wind. Therefore, it is interesting to consider 
whether the microphysical change in the strong field can affect the final results of the r- 
process nucleosynthesis. Furthermore, if we look into the scenario of hyperaccreting disks, 
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we should keep in mind that the disk will increase the wind entropy and affect its distribution 
in the polar region where propagates the stellar neutrino-driven wind (e.g., Nagataki et al. 
2006, 2007). Another possibility is that the disk outflows and the X-type winds, rather than 
the stellar neutrino-driven winds, will also produce the r-process elements with A > 130, 
because the neutrino-to-proton ratio in the disk and the outflows are sufficiently high, and 
the mass loss rate from the disk is much higher than that from the stellar surface (Kohri et 
al. 2005). However, whether the disk outflows can be the site for r-process nucleosynthesis is 
still under debate (Metzger et al. 2008a). Since the pressure of the outflow will drop quickly 
above the disk and be less than the magnetic pressure, the magnetic field above the disk 
may also play significant role in increasing the outflow entropy as discussed in Thompson 
(2003). As a result, more work should be done to understand the influence of a strong field 
on the r-process in the hyperaccreting magnetar systems. 



6.4. A Unified Scenario of Collapse-related-GRB Models 

Let us compare the hyperaccreting neutron star/mangetar model with the isolated mag- 
netar model, both of which are proposed to be able to produce GRB and GRB-like events. 
The hyperaccreting neutron star/magnetar disks can form in collapsars or compact binary 
mergers. In the collapsar scenario, rotation and magnetic fields make the core be possible 
to collapse into a massive neutron star/magnetar rather than a black hole. On the other 
hand, isolate magnetars are proposed to form via the rotating Type-Ib/c supernovae, the 
mergers of compact binaries, or the accretion-induced collapse of white dwarfs. Therefore, 
the progenitors of the hyperaccreting compact star model and the isolate magnetar model 
are actually very similar with each other. The main difference is the environment: the disk 
around the neutron star/magnetar will produce significant phenomena which will not be 
produced by the isolate magnetar. Let us see the case of massive star collapse for exam- 
ple. The isolate magnetar forms after a successful supernova explosion, and the supernova 
shock has created a cavity around the magnetar before it drives magnetically-dominated and 
Poynting-dominated winds. On the contrary, the unsuccessful or weak ongoing shock leads 
to the hyperaccreting system. The disk material around the neutron star/magnetar comes 
from the continuous infalling stellar envelope, or the fallback of stellar material which has 
been ejected by the shock but cannot reach sufficiently high velocity to escape the gravita- 
tional potential of the core. Kumar et al. (2008) discussed the accretion rate of different 
accretion stages, in which the accretion rate is different depending on which stellar zone is 
accreted. 

We try to use a unified point of view to consider the outcomes of the massive star 
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collapse. If the core collapse can initiate a successful supernova, then an isolated black hole 
or a rotating magnetar can form in the supernova remnant. Vietri & Stella (1998, 1999) 
discussed the possibility of a neutron star further losing its angular momentum by magnetic 
dipole or gravitational waves radiation during a long time and collapsing to an accreting 
black hole. In addition, the magnetar model shows the possibility that a variety of magnetic 
activities from the isolate magnetar can produce a GRB explosion after the supernova for 
a short time. On the other hand, if the outgoing shock formed during a core collapse 
cannot compete the continuous accretion from the stellar envelope, the collapse leads to the 
formation of a collapsar system, in which the type of the central compact object depends 
on many factors such as rotation, equations of state and so on. The black hole disk may 
produce a GRB explosion via neutrino annihilation or MHD processes, but the annihilation 
mechanism may not produce sufficient energy for energetic GRBs. The neutron star disk will 
increase the annihilation efficiency compared to its black hole counterpart, and the increased 
efficiency will be even higher if the central star is a magnetar. In the hyperaccreting magnetar 
system, the increased annihilation process and the magnetically-driven pulsar wind can work 
together to generate a more powerful jet than that generated by a single mechanism. The 
outflows and magnetically-dominated winds from the disk is possible to feed a late-disk- 
induced supernova associated with the GRB. In a word, the unified scenario shows that the 
outcomes have a closed relation to the initial stellar properties and the core collapse process 
itself. 

However, such a unified scenario as the GRB central engine candidate cannot be con- 
firmed directly. What we can observe is the GRB-related properties such as photon lightcurve, 
spectrum, neutrino emission, different bands of afterglows, which can be traced back to the 
central engine and show evidence for the existence of a central neutron star. We do not want 
to show further evidences here (see Dai 2004, Dai et al. 2006; Fan & Xu 2006; Yu & Dai 2007 
for more details). Finally we want to mention two points. One is that, the intermediate case, 
i.e., the fallback material forms a normal disk around a magnetar, has been discovered (Wang 
et al. 2006). Although the disk radiated IR emission is very different to the hyperaccreting 
disk, this discovery shows a link between isolate magnetar and hyperaccreting magnetar in 
the core collapses. Another point is that, since some stars that are less massive than Wolf- 
Rayet stars are more likely to form neutron stars and magnetars, we cannot preclude the 
possibility that the collapse of stars with an intermediate mass can also lead to the hyper- 
accreting neutron star/magnetar systems and produce the GRB-like phenomena. In fact, 
some evidences show that the GRB may be associated with Type II supernova (Germany et 
al., 2000; Rigon et al., 2003), although the evidences are still controversial today. 
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7. Conclusions 

The hyperaccreting neutron star or magnetar disk systems cooled via neutrino emission 
can form by the mergers of compact star binaries or the collapses of rotational massive stars. 
Strong fields of the magnetar can play a significant role in affecting the disk properties 
and even changing the accretion process. Our motivation in this paper is to investigate 
the influence of such strong magnetic fields on the disks, and discuss implications of the 
magnetar disk systems for the GRB and GRB-like events. Our conclusions are as follows. 

(i) We consider the magnetar field has a dipolar vertical component B z . The differential 
rotation between the disk and the magnetar will generate a toroidal field component B^, 
as well as a relatively weak radial component B r . The generated field can have an open 
or closed configuration, depending on the disk's viscous turbulence, magnetic diffusivity, 
disk angular velocity as well as twist limitation. Similar to pervious works, we use the 
parameter (3 to measure the strength of the toroidal field (equations (J3"Tj) and f[3"2l ). The 
generated large-scale disk field coupled with the accretion flows will transfer the angular 
momentum in radial direction and heat the disk via Joule dissipation together with viscous 
stress outside the Alfven radius r&- On the other hand, since the distribution and energy 
of electrons change significantly in the strong field environment, the disk pressure and a 
variety of neutrino cooling processes will be different compared to the case without fields. 
Therefore, we discuss the quantum effects of the strong fields on the disk thermodynamic and 
microphysical processes in Section 2, and list the MHD conservation equations to describe 
the behavior of the large-scale magnetic field coupling in the disk in Section 3. 

(ii) The quantum effects and field coupling in MHD equations play two competitive 
roles in changing the disk properties, the former to decrease the pressure, density and neu- 
trino luminosity with increasing field strength (Figure 4), while the latter to increase them 
(Figure 5). However, in most cases the large-scale field coupling is more significant than the 
microphysical quantum effect (Figure 7). Moreover, strong fields will change the electron 
fraction distribution Y e in the disk significantly. Larger peak of electron chemical potential 
i] e with more degeneracy electron states is obtained by stronger fields, while the change of 
r) e becomes more obvious for stronger fields. 

(iii) Similar to the neutron star disk without fields, the values of density, pressure, 
temperature, neutrino luminosity and electron chemical potential become higher for higher 
accretion rate, but the electron fraction Y e decreases with increasing the accretion rate. 

(iv) The magnetized disk maintaining a plane geometry would be more favorable for 
an open field configuration rather than a closed one. However, we still consider the two 
cases for completeness (Figures 6 and 9). For the disk with an open field and M = O.1M 
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s -1 , higher ratio of (3/s leads to higher density and pressure in the entire disk plane, higher 
temperature and neutrino luminosity in the inner part of the disk, as well as lower Y e in the 
outer part. Here s is the ratio of disk angular velocity and the Keplerian velocity. Also, 
electrons becomes more degeneracy at ~ 20 — 40 km for higher (3/s. On the other hand, the 
disk properties with a closed field not only depends on the values of (3 and s, but also the 
spin period of the central magnetar. Shorter period of the central star decreases the disk 
density, pressure, r] e , but increases Y e at the outer part of the disk. The effects of spin period 
are only obvious for sufficient field strength (e.g., ~ 10 16 G for the accretion rate O.1M s _1 ) 
from the central star. 

(v) The accretion flow in the disk plane outside the Alfven radius is viscously stable. 
However, whether the disk is thermally stable depends on many factors such as the disk 
region, magnetic field strength, disk angular velocity and accretion rate. Generally speaking, 
the disk will be definitely thermally unstable if its non-field disk counterpart with the same 
accretion rate is unstable. The disk region can also be thermally unstable even its non-field 
counterpart is stable, if the region is near the Alfven radius where magnetic field plays a 
more important role in transferring the angular momentum and heating the disk than the 
viscous stress. 

(vi) The thermally-driven outflows can also exist in the magnetar disk beyond the Alfven 
radius for the case of M < l.OM s _1 . We assume the accretion rate as a power law in radius 
for simplicity (Equation 045p ). The outflow will take away the disk angular momentum 
(Equation and may provide energy for a supernova associated with the GRB explosion. 
The disk density, pressure and neutrino luminosity decrease with increasing the outflow 
strength. Also, the ratio of magnetic and matter pressure -P;B/-Pmatter and the thickness 
of the disk become larger for stronger outflows. However, since the total energy taken by 
the thermal outflow depends on the size of the disk plane outside the Alfven radius, the 
total energy injection rate from the outflow decreases significantly for stronger stellar fields 
(Table [Q . Besides the thermally-driven outflow, strong fields inside the Alfven radius lead 
the accretion flow to co-rotating with the stellar field, which is possible to launch MHD 
winds along the field lines, and generate the X-type wind near the disk-magnetosphere 
boundary. These magnetically-dominated winds are nonrelativistic and may provide energy 
for a supernova explosion. 

(vii) In the hyperaccreting disks, the funnel accretion can only be important for ex- 
tremely strong fields, which depend on the accretion rate. The accretion process along the 
disk plane will be truncated in the stellar magnetosphere, so most of the accretion flow can 
be approximately considered as being lifted along the closed field lines onto the magnetar 
surface. In most cases, the hyperaccreting infalling funnel flow is unlikely to develop a shock 



-36 - 



with Mach number greater than unify (equation fl54l) ). The flow is accelerated by the grav- 
itational force and transfers the gravitational binding energy to the kinematic energy and 
next to the heating energy near the magnetar surface. The funnel flow will cover a ring-like 
belt of "hot spot" around the magnetar surface and emit thermal neutrinos. Because the 
temperature is higher and neutrino emission region is more concentrated than the hyper- 
accreting neutron star without fields, the funnel accretion can accumulate more powerful 
neutrino annihilation luminosity (Figure 13). 

(viii) The neutrino annihilation process both from the magnetar surface and from the 
disk plane will be higher than that without fields. Moreover, the neutrino annihilation mech- 
anism and the magnetic activity from the stellar surface (i.e., the pulsar wind mechanism) 
can work together to generate and feed an ultra-relativistic jet along the stellar magnetic 
poles. If the stellar spin period is sufficiently short (e.g., ~ 4 ms for the field ~ 10 16 G and 
M = O.1M s" 1 ), the jet from the magnetar will be magnetically-dominated and mainly 
feeded by extraction the stellar rotational energy. If the magnetar spin period is long (i.e., 
longer than the critical value P rc in Equation (|62j) ). the jet is thermally-driven and feeded 
by the annihilation process. In the intermediate case, on the other hand, the relativistic jet 
can be launched by the pulsar-wind-like process and neutrino annihilation together. 
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A. Disk Equations without Field 

The thermodynamical, microphysical and conservation equations in a hyperaccreting 
disk without fields can be seen in many previous works. In order to compare them with the 
case with strong fields, we list them systematically in this appendix. Here we do not consider 
equations of the neutron star inner disk with a self-similar structure as discussed in paper I 
and paper II. 
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The electron/positron density number reads 



_ (m e c) 3 r°° x 2 dx 
and the electron/positron pressure is 

?± _im^ r^; dx 

e± 3n 2 h 3 J ^/x 2 + 1 e (mec 2 v^+T TMe )/fcBT + ! ' 1 J 

The total pressure is the summation of the pressure of electrons, nucleons, radiation and 
neutrinos (without magnetic pressure): 

P = P e - + Pe+ + Pnuc + Pad + P„, (A3) 

where P nuc = pRT and P ra( j = asT 4 /3 with i? being the gas constant and as being the 
radiation constant. 

The total neutrino cooling rate is the same as equation ([9]) with various absorption and 
scattering depths showed in the beginning of Section 2.2. 

The neutrino cooling by electron-positron captures by nucleons are contributed by three 
terms as follows: 

/■oo 

q p+e -^ n+Ve = Kn p / e{e 2 - l) 1/2 {e - q) 3 f e ~de, (A4) 

Jq 

/oo 
e{e 2 - l) 1/2 (e + q) 3 f e +de, (A5) 

gn^ P+e -+p e = #n n ^\(^ 2 - l) 1/2 (q - e) 3 (l - / e _)de, (A6) 

when the field is absent or can be ignored. Here K, e, q and / eT as in Equations (ITTj) to (TT5|) 
in Section 2.2. Other cooling rates without fields have been showed in Section 2.2. 

The chemical equilibrium is 

ln (-) = /(r.)^r^ + 1 1 - /^)]^r^. ( A7 ) 



rip/ &bT fcgT 



where the weight factor f(r u ) = exp(— t v ) combines the formula from the neutrino-transparent 
limit case with the neutrino-opaque limit case of the /3-equilibrium distribution, and Q = 



qm e c 2 . 



A set of conservation equations without fields are as follows: 
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Mass conservation (continuity) equation reads 

M = -27rrSw r , (A8) 

this equation keeps the same with and without fields in the vertically-integrated disks. How- 
ever, we should consider another continuity equation for outflows and winds as 

M = -47rr 2 pv r (An/n), (A9) 

where Afl is the opening solid angle of the wind. 

The integrated angular momentum conservation equation is 

M 

3^/ = i/E, (A10) 
The local energy conservation without fields is 

Qt* = Q^ + QZ, (Ail) 

where the viscosity heating term is 

Q^^f (A12) 

and the advection term is 



R M _ 4 aT 3 
- (1 + Y e ) + -g. 



(A13) 



2 v w y ~ p 

with g* — 2 for photons and g* — 11/2 for a plasma of photons and relativistic e~e + pairs 
The charge conservation equation is 

pYe 



= n p = n e - - n e +, (A14) 
m B 

which will not be changed in strong fields if we do not consider other charged particles 
created in the strong field environment. 



B. Solid Angle of Funnel Flow 

We take two methods to estimate the latitude and area of the ring-like belt of the "hot 
spot" on the magnetar surface. The first one is just to follow the approximation made in 
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Section 5.2. We assume that the funnel accretion flow is lifted in the region with radius 
between r m and r^, and the disk plane will be truncated at the magnetosphere r m . Taking 
r m — r A — 8, we have the field line equation from the two disk radius r A and r m as 

rn = (Bl) 
sin 6 



ra = (r A -5)^, (B2) 
sin B 



or we obtain 



sinfl, = ( — ) sinB, (B3) 



rA 



1/2 



sinfl 2 = sinf^i + AO) = =— sinB, (B4) 

If 5 -C rA, sinB ~ 1, we can obtain an analytic solution of AO and AQ/2tc. We have 

sin# 2 = sinfli + cos^Atf = (^j (l + , (B5) 

or 

M = tt4 V (B6) 

1 \/r A \r A - n) 

The solid angle spanned by the hot spot ring-like belt can be calculated as 

p2ir p62 

AQ= d<p sinOdO (B7) 
Jo Je 1 

and we obtain an analytic result of the ratio Af2/27r as a function of r A and 5. 

AQ 27rsin0iA0 5 I F* , , 

(B8) 



27r 27T 2r A y r A -r* 

We do not want to solve detailed conservation equations or a set of MHD differential equa- 
tions as in Ustyugova et al. (1999). Table [3] gives the results with different parameters of 
r A /r*, and B as in the case of hyper accreting disks with strong fields. The typical value 
for An/27r is around 10~ 2 , except for the case when r A is small and 5 ~ r A , which can still be 
considered as disk plane accretion. Stronger field makes the funnel flow be more significant 
and the value Af2/27r be smaller. 

Next we give another scenario based on the consideration that the disk has a thickness 
near the magnetosphere r m . We try to give a simple mathematical model: The ring-like belt 
is formed by the accretion matter at r m with different height z in the vertical direction. In 
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other words, the ring-like belt area on the magnetar surface can be traced back to the disk 
plane vertical section at r m . 

sm 2 9 

ni = r m —^-, (B9) 
sin 6 

, sin 2 ^ 

ri2 = r'-^-, (BIO) 
sin B 



where r' satisfies 



cos 2 (if/r m ) 

Similarly, if sinB ~ 1 and H 2 /r 2 ^ <C 1, we have 

1/2 



sin 2 6, (Bll) 



sin^ = ( — ) sinB, (B12) 



A^^/^f-V. (B14) 

2 



4 V r m — r* V r 



AQ _ 27rsingi A9 _3 [77 I K f H_ \ (m5) 
2vr 27T 4y r m \j r m - r* \r m J 

We can take fraction of r^. Table S] gives numerical results in this case. The typical 

value of AQ/2tt is as ~ 10 -3 — 10~ 2 , which can be slightly smaller than that obtained using 
the first estimation. 

On the other hand, the accretion rate onto the neutron star without fields, as discussed 
in paper II, can be estimated as (27r7\,,if )/(27rr 2 ) = H/r* (note that we consider the disk 
over a half thickness). Thus the solid angle in this case is in the order of > 0.1, which is 
larger than the area formed by funnel accretion. 
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Disk Heating Rate and Maximum Outflow Energy Rate 







Disk Heating Energy 






Max Outflow Energy 








Rate (10 51 erg s" 1 ) 






Rate (10 51 erg s" 1 ) 






B = 10 15 G 


Bo = 10 16 G 


Bo = 10 17 G 


B Q = 10 15 G 


Bo = 10 16 G 


Bo = 10 17 G 


£=0.2 


29.3 


25.9 


8.08 


8.60 


8.63 


1.45 


£=0.6 


18.3 


15.7 


6.47 


19.5 


18.8 


3.07 


£=0.9 


13.5 


11.7 


5.53 


24.4 


22.9 


4.01 



Table 1: The disk heating rate is calculated in the region outside the Alfven radius, where 
differential rotation and viscosity are significant: -Eheat — X-^" 1 ' 27rr<ir. Here the outflow index 
£ = 0.2,0.6,0.9, and the accretion rate M = O.5M s _1 . The max thermally-driven energy 
rate is estimated using equation f|4T|) . The total thermal outflow energy can be considered 
as 0.1 — 1 fraction of the maximum energy, as discussed in paper II. 



Disk Properties Depending on Magnetar Field Strength 



B - 


- 10 14 


G 


no funnel accretion, disk is similar to that without fields 


B - 


- 10 15 


G 


no funnel accretion, MHD coupling is important in the disk inner region, thermal outflow 


B - 


- 10 16 


G 


weak funnel accretion, disk is denser, hotter with higher pressure, brighter L„ p from the disk, 
thermal outflows from r > Va, magnetic winds inside 


B - 


- 10 17 


G 


strong funnel accretion, but no shock, much brighter L„ p from the stellar "hot spot" , 
magnetically-dominated wind significantly 



Table 2: The accretion rate near the Alfven radius ta or the stellar surface r* (depending on 
ta > r* or < r* respectively) is M = O.1M s -1 . The neutrino cooling emission is efficient 
in this case, and the accretion flow is an NDAF. The inner region satisfying the self-similar 
structure (papers I and II) can be ignored even the disk is similar to that without fields. 
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Values of Af2/27r on Magnetar Surface Based on the First Scenario 



5/n 


2 


r m /r*(9 = 90°) 
4 


7 


2 


r m /n (9 = 75°) 
4 


7 


0.1 


1.84e-2 


6.38e-3 


2.74e-3 


1.77e-2 


6.16e-3 


2.64e-3 


0.5 


0.110 


3.46e-2 


1.43e-2 


0.106 


3.34e-2 


1.38e-2 


1 


0.292 


7.74e-2 


3.02e-2 


0.282 


7.48e-2 


2.92e-2 


2 




0.208 


0.070 




0.200 


6.70e-2 



Table 3: The ratio of the solid angle of the "hot spot" ring-like belt AQ formed by funnel 
accretion to the half total solid angle 2tt based on the first scenario in Append B, in which 
the accreted matter onto the ring-like belt of "hot spot" is from the region between the 
magnetosphere edge r m and the Alfven radius ta- We take from 0.1 to 2, and r m /r* 
from 2 to 7. Here the left three columns are for 9 = 90° and the right three ones are for 
9 = 75°. The typical value of AQ is around 10~ 2 . 



Values of Af2/27r on Magnetar Surface Based on the Second Scenario 



H/r m 


2 


9 = 90° {r A /n) 
4 


7 


2 


9 = 75° {r A / n ) 
4 


7 


0.1 


5.28e-3 


3.74e-3 


2.82e-3 


5.50e-3 


3.62e-3 


2.72e-3 


0.3 


4.60e-2 


3.26e-2 


2.46e-2 


4.44e-2 


3.14e-2 


2.38e-2 


0.5 


0.120 


8.50e-2 


6.42e-2 


0.116 


8.22e-2 


6.20e-2 



Table 4: The ratio of the solid angle of the "hot spot" ring-like belt AQ to the total solid 
angle 2n based on the second scenario in Append B, in which the accreted matter onto the 
ring-like belt of "hot spot" is from the disk vertical section surface at r = r m . We take the 
disk thickness at this radius as H/r m = 0.1, 0.3, 0.5, and r m /r* also from 2 to 7. The typical 
value of Ail is still around 10~ 3 — 10~ 2 in this scenario. 
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Fig. 1. — Alfven radius function of accretion rate in hyperaccreting disks. We take 

the magnetar field as 10 15 , 10 16 and 10 17 G. 



-50 - 




Fig. 2. — Left: Convergence performance of Landau level series of electron density and 
pressure with different magnetic fields B = 10 16 G (solid line), 10 15 G (dashed line), 10 14 G 
(dotted line) and 10 13 G (dash-dotted line) and temperature T = 5 x 10 10 K and chemical 
potential rj e = 1, where b n is the n term in the density and pressure series © and (J5J). Right: 
Comparison of density and pressure with and without magnetic fields with temperature 
T = 10 10 K (solid line), 5 x 10 10 K (dashed line) and 2 x 10 11 K (dotted line). 
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Fig. 3. — Left: Convergence performance of Landau level series of neutrino cooling with 
different magnetic fields. Upper one for the cooling rate q p + e -^n+v e and lower for g n+e +_> p+ p e , 
where temperature T = 5 x 10 10 K and chemical potential r] e = 1 is adopted. Lines are as 
in Fig. 1. Right: Comparison of neutrino cooling rates with and without magnetic fields. 
Temperature T = 5 x 10 10 K and rj e = 0.1, 1 and 10. Also, upper is for q p+e -^ n+Ue and lower 

for Qn+e+^p+Pe- 
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Fig. 4. — Quantum effects of microphysics and thermodynamics in strong magnetic fields 
to affect the disk properties. In order to see the effects clearly, we still take the non- 
magnetized conservation equations (i.e, ( 1A8j) . ( 1A10I) to ( 1A14I) in Appendix A), and only 
change a set of equations of state in Section 2. Also, we take the magnetic field as uniform 
B = 10 14 , 10 15 , 10 16 , 10 17 G and without field, the accretion rate M = O.1M s _1 . 



-53 - 




Fig. 5. — Effects of strong magnetic field coupling on angular momentum transfer and heating 
via Joule dissipation on the disk. We do not consider the field quantum effects as discussed 
in Section 2, but only adopt the MHD equations in Section 3. The magnetic field topology 
is dipolar, and an open configuration with parameter (3 = 0.6, angular velocity ratio s — 1, 
and B = 1 14 , 10 15 , 10 16 , 10 17 G. The accretion rate is the same as in Figure 4. 




Fig. 6. — Disk structure with M = O.1M s^ 1 for the magnetar surface vertical field B = 
10 16 G. The magnetic field in the disk is open with parameter (3 = 0.2, 0.6, 1 and disk angular 
velocity s = 0.5, 1. 
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Fig. 7. — Disk structure with M = O.1M s^ 1 for the magnetar surface vertical field B = 
10 14 , 10 15 , 10 16 and 10 17 G, where the magnetic field in the disk is in an open configuration 
with (3 = 1 and s — 1. 
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Fig. 8. — Disk structure for different disk accretion rate M = 0.02 (solid line), 0.1 (dashed 
line) and 1 M s" 1 (dotted line) with open disk field (3 = 0.6 and s — 1, while the surface 
vertical field B = fO 16 G. 
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Fig. 9.— Disk structure with M = 0.1M© s" 1 for the surface vertical field B = 10 15 , 10 16 
G, the disk magnetic field is closed with j3 ~ 1, s — 1. The period of a central magnetar 
P r = 0.005,0.01 and 0.1 sec. 
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Fig. 10. — M — E curves for given disk radius r =40 and 80 km. In the case of an open disk 
field we take f3 = 0.6 and s — 1; while in the case of a closed disk field we take /3 — 1, s — 1 
and P r = 100ms. 
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Fig. 11. — Ratio of magnetic pressure P B to the accretion matter P ma tter, and ratio of 
neutrino emission Q~ to the local heating rate Q^ is + <5^ ul as functions of disk radius. 
Upper two panels for the accretion rate M = O.1M s" 1 . Open disk field with B = 10 16 
G, j3 — 0.2, s = 1 (thick solid line), j3 = 0.6, s = 1 (thick dashed line), j3 — l,s — 1 
(thick dotted line), B = 10 17 G, f3 = l,s = 1 (thick dot-dashed line) and closed field with 
j3 = l,s = l,P r = 0.005 s (thin solid line). Bottom two panels for the accretion rate 
M Q = 1M Q s^ 1 , and the lines are the same as the upper panels. 
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Fig. 12. — Disk structure and luminosity for three values of outflow index £ = 0.2, 0.6, 0.9 
with the open disk field (3 = 0.6, s — 1, B — 10 16 , 10 17 G. We give an initial accretion rate 
M = 0.5Mp, s _1 at the outer radius r = 124 km. 
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Fig. 13. — The neutrino annihilation rate l vV (erg s 1 cm 3 ) along the stellar mag- 
netic pole as a function of height, where the values of sin# ~ ^/t*/ta are taken as 
sin6 = 0.2,0.4,0.6,0.8,1.0, Afi/27r = 10~ 2 , accretion rate M = O.5M s _1 , and the 
energy release efficiency e = 0.5. The integrated annihilation energy along the pole is 
J^Udz = 4.95 x 10 41 ,3.31 x 10 40 ,5.92 x 10 39 , 1.31 x 10 38 and 9.21 x 10 37 erg s' 1 cm' 2 
for the value of sin^ = 0.2, 0.4, 0.6, 0.8 and 1.0 respectively. 



-62 - 




Fig. 14. — Neutrino luminosity for accretion rate M = 0.01, 0.1 and 1 M s _1 and magnetar 
surface field B = 10 15 , 10 16 and 10 17 G without Joule dissipation in the disk. The open field 
model is adopted with f3 = 0.6 and s — 1. 



